arXiv:1502.06138vl [math.SP] 21 Feb 2015 


Rational invariant tori and band edge spectra for 
non-selfadjoint operators 


Michael Hitrik 
Department of Mathematics 
University of California 
Los Angeles 
CA 90095-1555, USA 


Johannes Sjostrand 
IMB, Universite de Bourgogne 
9 avenue Alain Savary, BP 47870 
21078 Dijon cedex, France 
and UMR 5584 CNRS 

j ohannes. sj ost r and@u-bour gogne. fr 


hitrik@math.ucla.edu 


Abstract 


We study semiclassical asymptotics for spectra of non-selfadjoint pertur¬ 
bations of selfadjoint analytic h-pseudodifferential operators in dimension 2, 
assuming that the classical flow of the unperturbed part is completely in- 
tegrable. Complete asymptotic expansions are established for all individual 
eigenvalues in suitable regions of the complex spectral plane, near the edges 
of the spectral band, coming from rational flow-invariant Lagrangian tori. 
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1 Introduction 

Spectra for semiclassical non-selfadjoint operators often display fascinating fea¬ 
tures, from lattices of low-lying eigenvalues for operators of Kramers-Fokker-Planck 
type a. 0 to eigenvalues for operators with analytic coefficients in dimension one, 
concentrated to unions of curves, Ell. Ba. 0.0. The work |20] has established 
that for wide and stable classes of non-selfadjoint analytic pseudodifferential op¬ 
erators in dimension two, the individual eigenvalues can be determined accurately 
in the semiclassical limit, by means of a complex Bohr-Sommerfeld quantization 
condition, and form a distorted two-dimensional lattice. Now in many natural situ¬ 
ations [H], [2^, [26], [28], one encounters non-selfadjoint operators of the form 

Pe=p{x,hD^) + ieq{x,hD^), 0 < £ < 1 , ( 1 . 1 ) 

considered on or a compact real analytic manifold, with Pe=Q being selfadjoint. 
Here 0 < h 1 is the semiclassical parameter and the second small parameter e 
represents the strength of the non-selfadjoint perturbation. The principal symbol 
of Pg in fll.ip is of the form Pe{x,^) = p{x,^) -1- ieq{x,^), where p is real, and let 
us also assume, to fix the ideas, that q is real. Both p and q are assumed to be 
analytic, with p elliptic near infinity. The spectrum of Pg near the origin is confined 
to a band of width 0{e), and the general problem is to understand the distribution 
of eigenvalues of Pg near 0, in the semiclassical limit h —>■ O’*". To this end, let us 
assume that 0 is a regular value of p, so that the energy surface p“^(0) is a smooth 
compact submanifold of the phase space. We then know [m, ng that the real parts 
of the eigenvalues of Pg near 0 are distributed according to the same Weyl law as 
that for the unperturbed operator Pe=o. In order to study the distribution of the 
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imaginary parts of the eigenvalues, following the method of averaging Ell. IB , we 
let Hp be the Hamilton vector held of p and introduce the time averages 


1 

{q)T = 7 p g o exp (tHp) dt, T > 0, (1.2) 

J J-T/2 


of q along the ifp-trajectories. It follows from na, IZ71. ca that if 2 ; G Spec(Pe) is 
such that |Re 2 ;| < 6, then 


Im z 

lim inf {q)rp — o{l)< - < lim sup + o(l), (1.3) 

T^oop-i(O) e T^oop-i^ 


as (e, 6, h) —)■ 0 +. 

The spectral analysis for non-selfadjoint operators of the form (11.11) has been pur¬ 
sued by the authors in the series of papers 0- iia. the latter work jointly with 
S. Vu Ngoc, when the dimension n = 2 and the Hp-fiow is either periodic or com¬ 
pletely integrable. Let us focus, from now on, on the completely integrable case, 
which will be considered also in the present work. In this case, the energy surface 
p“^(0) is foliated by invariant Lagrangian tori, along with possibly some other more 
complicated flow-invariant sets. When A C p“^(0) is an invariant torus such that 
the rotation number of Hp along A is Diophantine, i.e. poorly approximated by 
rational numbers, or more generally, irrational, we have that the time averages {q)rp 
along A converge to the space average (g)(A) of q over A, as T ^ 00 . When A is 
a torus with a rational rotation number, or a singular set in the foliation of p~^( 0 ), 
then in analogy with (11.31) . we introduce the compact interval 


Qoo(A) 


lim inf(g).^, lim sup{q)rp 

T—>00 A T—)-oo ^ 


(1.4) 


of limits of the time averages above. 

Let Fq ^ R- be such that Fq = (g)(Ad) for a single Diophantine Lagrangian torus 
Ad C p“^(0), and let us assume that 


Fq ^ Qoo(A), (1-5) 

for any other invariant set A 7 ^ A^ in p“^(0). It was then shown in |T2] that 
the spectrum of can be determined completely, modulo 0{h°°), in a rectangle 
of the form [—h^/C, h^/C] -|- ie[FQ — /C, Fq + /C], where 5 > 0 and £ satishes 

< e ^ 1, for K ^ 1. Similarly to [20], the spectrum has a structure of a distorted 


3 








two-dimensional lattice, with the horizontal spacing ~ h and the vertical one ~ eh. 
A closely related result was obtained in m, giving a Weyl type asymptotic formula 
for the number of eigenvalues of in an intermediate spectral band, bounded from 
above and from below by Diophantine levels, such as Fq above. It turned out that 
the distribution of the imaginary parts of the eigenvalues of is governed by a 
Weyl law, expressed in terms of phase space volumes associated to p and the long 
time averages of q. 

Having elucidated the role played by flow-invariant Diophantine tori in the spectral 
analysis of Pg, let us now turn the attention to spectral contributions of tori that 
are rational, which constitutes the subject of the present work. Let Pq G R be such 
that Pq = {q){Ad), for a Diophantine torus as above, and rather than demand¬ 
ing (II.5j) . let us assume that there exists a rational torus A^ C p“^(0) such that 
Po e Qoo(Ar), Po 7 ^ (<?)(Ar), while Pq ^ Qoo{A), for A ^ Ad, A^. An attempt to 
determine the individual eigenvalues of Pg near iePo was made by the authors in the 
work [To], by means of the normal form techniques. As a result, the normal forms 
near A^ that we obtained were given by a family of one-dimensional ’’resonant” non- 
selfadjoint operators, and the possibility of quite serious pseudospectral phenomena 
for this family [2] prevented us from computing the eigenvalues individually. Cor¬ 
respondingly, the main result of jTOl was weaker, establishing that the spectrum of 
Pg near iePo was of the form EdU Er, where the ’’Diophantine” contribution Ed is 
a distorted lattice that can be described explicitly, as in [12] , and the cardinality of 
the ’’rational” contribution Er is <C than that of Ed- 

Subsequently, in the course of some numerical experiments, the authors have en¬ 
countered peculiar pictures of the spectra of Pg, where the eigenvalues had the form 
of a ’’centipede”, with the body agreeing with the range of torus averages of g — 
see Section [H] for the illustrations and the details of the numerical computations. 
The legs of the centipede were more mysterious at hrst, but things became clearer 
when we realized that they represented the influence of suitable rational tori. It 
became then natural to hope that the eigenvalues near the extremities of the legs 
could be determined asymptotically in a rigorous way, since the pseudospectral ef¬ 
fects should become more moderate near the edges of the spectral band, [2]. The 
main result of the present work, giving a complete asymptotic description of the 
individual extremal eigenvalues of Pg, can be considered as a justihcation of this 
hope. 

Let us conclude the introduction by formulating, in a rough way, the main result of 
the paper — see Theorem 2.1 below for the precise statement. Let Aq C p“^(0) be 
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a rational Lagrangian torus such that 

infQoo(Ao) < inf (infQoo(A)). (1.6) 

A^Ao 

The restriction of the ifp-flow to Aq is periodic with primitive period Tq > 0, and 
the time average {q)rp^ in (ll.2p can naturally be viewed as a function on the space 
of closed orbits Ao/exp(Rihp). Let us assume that viewed as a function on 

Ao/exp(RiLp), has a unique minimum which is nondegenerate, and restrict £ to a 
suitable interval of the form rj > 0. Then for any hxed Co > 0 the 

eigenvalues of in the region 

I h Tm 2^ hi 

\z eC; |Re z\ < , -< infQoo(Ao) + Cq—^ \ 

[ Co^/e e y/e) 

can be determined completely, modulo 0{h°°), and are given by 

1 iS* 

Ayfc = a(^2) + ieb{^2) + 6 = h{j “ j eZ, keN. (1.7) 

Here a(0) = 0, a'(0) > 0, 6(0) = inf(5oo(Ao), and S and ko are the classical action 
and the Maslov index of a primitive closed Rp-trajectory in Aq. We have a complete 
asymptotic expansion for ^ in integer powers of 6. = 

OO 

u=0 

where 

A^(0,0) = de*"/^(2fc + l), d>0. 

Acknowledgement. We are very grateful to Michael Hall and San Vu Ngoc for 
interesting discussions concerning numerical computations of eigenvalues and for 
showing us promising numerical results for a damped spherical oscillator. The second 
author is supported by the project NOSEVOL ANR 2011 BS 01019 01. 

2 Statement of the main results 

2.1 General assumptions 

We shall start by describing the general assumptions on our operators, which will 
be the same as in [TO] , na, as well as in the earlier papers mentioned above. Let 
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M denote either the space or a real analytic compact manifold of dimension 2. 
When M = R2, let 

Ps = P^{x,hDx,e;h), 0 < h < 1, (2.1) 

be the h-Weyl quantization on R^ of a symbol P{x, e; h) (i.e. the Weyl quantiza¬ 
tion of P{x, h^, E] h)), depending smoothly on £ G neigh(0, R) and taking values in 
the space of holomorphic functions of (x, in a tubular neighborhood of R^ in C^, 
with 

\P{x,^,£;h)\ < C>(l)m(Re(x,0), (2.2) 

there. Here 1 < m G (^““(R^) is an order function, in the sense that 

m{X) < Co{X - Y)^°m{Y), X,Y e R^ (2.3) 

for some Cq, Nq > 0. We shall assume, as we may, that m belongs to its own symbol 
class, so that = Oaim) for each a G N^. Then for h > 0 small enough and 
when equipped with the domain H{m) := {mP{x^hD))~^ (L^(R^)), P^ becomes a 
closed densely dehned operator on L^(R^). 

Remark. The analyticity assumptions will allow us to treat the case when e: x 
for 0 < 5 < 1. When e: = 0{h), standard C'°°-microlocal analysis would have been 
sufficient. 

Assume furthermore that 


F(x, e;h) 0 (2-4) 

j=0 

in the space of holomorphic functions depending smoothly on e: G neigh(0, R) and 
satisfying fl2.2p in a fixed tubular neighborhood of R"^. We assume that po,e is elliptic 
near inhnity, 

\Po,sAA)\> ^m{Re{x,A), |(x,OI>C', (2.5) 

for some C > 0. 

When M is a compact manifold, for simplicity we shall take P^ to be a differential 
operator on M, such that for every choice of local coordinates, centered at some 
point of M, it takes the form 

= X] aa,£(a;;^)(^^x)", (2.6) 

\a\<m 
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where tta^eix; h) is a smooth function of e G neigh(0, R) with values in the space of 
bounded holomorphic functions in a complex neighborhood of a; = 0, independent 
of h when |q;| = m. We further assume that 

OO 

~ X] <^a,e,j{x)h\ /l 0, (2.7) 

i=o 

in the space of such functions. The semiclassical principal symbol po,£) dehned on 
T*M, takes the form 

Po,s{^^ 0 = X] aa,e,o{^)C, (2.8) 

if {x,^) are the canonical coordinates on T*M. We make the ellipticity assumption, 

ipo,=n.oi>^«r. {x,()erM, i«i>c, (2.9) 

for some large C > 0. Here we assume that M has been equipped with some real 
analytic Riemannian metric so that |.^| and {^) = (1 + are well-dehned. 

Sometimes, we write for po,e and simply p for po.o- We make the assumption that 

P£=o is formally selfadjoint. 

In the case when M is compact, we let the underlying Hilbert space be fi{dx)) 

where fi{dx) is the Riemannian volume element. 

The assumptions above imply that the spectrum of Pg in a hxed neighborhood of 
0 G C is discrete, when 0 < h < ho, 0 < e < Eq, with ho > 0, Sq > 0 sufficiently 
small. Moreover, ii z & neigh(0, C) is an eigenvalue of then Im;^ = 0(e). 

We furthermore assume that the real energy surface p~^(0) r\T*M is connected and 
that 

dp 0 along p“^(0) fl T*M. 

In what follows we shall write 


Pe = p + ieq + 0{e^), (2.10) 

in a neighborhood of p~^(0) fl T*M, and for simplicity we shall assume throughout 
the paper that q is real valued on the real domain. (In the general case, we should 
simply replace q below by Re q.) We let Pdp = p'^ ■ dx — p'^ ■ be the Hamilton vector 
held of p. 
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2.2 Assumptions related to the complete integrability 

As in na, |in] . let US assume that there exists an analytic real valued function 
/ near p“^(0) fl T*M such that Hpf = 0, with the differentials df and dp being 
linearly independent on an open and dense set C neigh(p“^(0) r\T*M,T*M). For 
each E G neigh(0, R), the level sets Aa^E = /~^(a) r\p~^{E) fl T*M are invariant 
under the ifp-flow and form a singular foliation of the 3-dimensional hypersurface 
p~^{E) n T*M. At each regular point (i.e. non-critical point for the restriction 
of / to p~^{E)), the leaves of this foliation are 2-dimensional analytic Lagrangian 
submanifolds, and each regular leaf is a hnite union of tori. In what follows we shall 
use the word “leaf” and notation A for a connected component of some Aa^E- Let J 
be the set of all leaves in p“^(0)nT*M. Then we have a disjoint union decomposition 

p-^(O) nT*M = y A, (2.11) 

AeJ 

where A are compact connected Rp-invariant sets. The set J has a natural structure 
of a graph whose edges correspond to families of regular leaves and the set S of 
vertices is composed of singular leaves. The union of edges J\S possesses a natural 
real analytic structure and the corresponding tori depend analytically on A G J\S 
with respect to that structure. See section 7 in [10] for an explicit description of the 
Lagrangian foliation in the case when M is an analytic surface of revolution in R^. 

In what follows, we shall assume that the graph J is hnite. We shall identify each 
edge of J analytically with a real bounded interval and this determines a distance 
on J in the natural way. Assume that the following continuity property holds. 

For every Aq G J and every £ > 0, there exists 5 > 0, such that if (2.12) 
A G J, distj(A, Aq) < 5, then A C {p G p~^(0) fl T*M; dist(p, Aq) < e}. 

Remark. Let us assume that / is a Morse-Bott function when restricted to p~^(0) fl 
T*M, in the sense that the set of critical points of the restriction of / to p“^(0)nT*M 
is a disjoint union of connected submanifolds, with the transversal Hessian of / 
being nondegenerate along each of the submanifolds. In this case, the structure of 
the singular leaves is known [30|. The set J is then a hnite connected graph and the 
property fl2.12p holds. 

Each torus A G J\S carries real analytic coordinates Xi, X 2 , identifying A with 
T2 = R727rZ2, so that along A, we have 

Rp = + 02^X2, (2.13) 
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where Oi, 02 G R. The rotation number is dehned as the ratio 


a;(A) = [oi : 02 ] G RP^, 


and it depends analytically on A G J\S. Recall also that the leading perturbation q 
has been introduced in fl2.10l) . For each torus A G J\S, we dehne the torus average 
(g)(A) obtained by integrating q\\ with respect to the natural smooth measure on 

A. 


We introduce the time averages, 

1 

{q)T= 7 ^ qoexp{tHp)dt, T > 0, (2.14) 

I J-T/2 

and consider the compact intervals Qoo(A) C R, A G J, dehned as in [T2] . 


Qoo(A) 


lim inf(g)r, lim sup(g)r 

T—^■oo A T—^00 yY 


( 2 . 15 ) 


Notice that when A G J\S and ^(A) ^ Q then Qoo(A) = {(g)(A)}. In the rational 
case, we write Ci;(A) = where m G Z and n G N are relatively prime, and where 
we may assume that m = 0{n). When k{uj{A)) := \m\ + \n\ is the height of Ci;(A), 
we recall from Proposition 7.1 in [12] that 


Q.(A)c(,)(A) + o(jj^) [-1,1], 


(2.16) 


Remark. As J\S 9 A —)■ Aq G S', the set of all accumulation points of {q){A) is 
contained in the interval Qoo(Ao). See the related remark in [10], Section 2. 

From Theorem 7.6 in [12] we recall that 


1 

£ 


Im {Spec{Pe) n {z] |Re^| < 5}) C 


inf IJ Qoo(A) 
AeJ 


as (e, h, 5) —)■ 0. 


o(l),sup IJ Qoo(A) + 0 ( 1 ) , 
AeJ 

(2.17) 


2.3 The main result 

Let Af) e J\S be a rational invariant Lagrangian torus, so that as above, too 
ttlAo) = 5 € Q, m = 0{n). Assume that the isoenergetic condition holds, 

(dAa;)(Ao) ^0. (2.18) 
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We recall from Section 2 of [TO] the behavior of the interval Qoo(A) when A 7 ^ Aq is 
a rational torus in a neighborhood of Aq. Writing a;(A) = | where p G Z and g G N 
are relatively prime, p = 0{q), we get, using that a;(A) 7 ^ wq, 


a;(A) — cuol > — > 
nq 


1 

nfc(a;(A)) ’ 


(2.19) 


and therefore, in view of (I2.16I1 . 

Qoo(A) C (g)(A) + 0(dist(n;(A),a;o)°°)[-l,l]. (2.20) 

This estimate is uniform in uq provided that we have a uniform upper bound on the 
height of the rotation number oiq G Q. 

Let us assume that the chosen rational torus Aq is such that 

infQoo(Ao)< inf infQoo(A). (2-21) 

AeJ\{Ao} 

The result below remains valid with the obvious modihcations, if we replace fl2.2ip 
by the assumption 

supQoo(Ao) > sup supQoo(A). (2.22) 

AeJ\{Ao} 

Let us choose, as we may, action-angle coordinates (a:, 0 Aq, so that Aq is given 
by = 0} in x R|, p = p{^), and so that 

p{0) = 0, %p(0) = 0, %p(0) > 0, dlp{0) ^ 0. (2.23) 

Here the last property follows from fl2.18l) . and in order to £x the ideas, we shall 
assume that 9|^p(0) > 0. By the implicit function theorem we have 

%p(O = 0^6 = /(6), (2.24) 

where / is an analytic function with /(O) = 0, and we obtain an analytic family of 
rational Lagrangian tori A^ C p~^{E), E G neigh(0, R), given by 

^2=UE), = (2.25) 

Here ^2 = ^ 2 {E) is the unique smooth solution of the equation p{f{^ 2 ),^ 2 ) = E, 
close to 0, such that ^2(0) = 0. 
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Writing q = q{x,^) in terms of the action-angle coordinates (x,^), let 


(?) 2 (a;i ,0 = ^ / (lix: 0 dx 2 , ^ e neigh( 0 , R^), (2.26) 


1 


be the average of q with respect to X 2 . We assume that 


T 9 xi {q) 2 {xi, 0) has a unique minimum 
which is nondegenerate. 


(2.27) 


In order to give an invariant description of the assumption 02.271) , let us notice that 
when restricted to Aq, the Hamilton flow of p is periodic of primitive period Tq > 0 
and the average (g) 2 (xi, 0 ) can naturally be viewed as the flow average {q)rp^, dehned 


as in 02.14p . considered as a function on the space of closed iTp-orbits in Aq, 

Ao/exp {KHp) ~ T. 


In its invariant formulation, the assumption 02.27P therefore states that flow aver- 
age {q)rp^, viewed as a function on Ag/exp (RRp), has a unique minimum which is 
nondegenerate. 

It follows from 02.271) that the function T 3 Xi i—)■ ( 5 ) 2 (xi, ^) has a unique minimum 
Xi = Xi(^) which is nondegenerate, for ^ G neigh(0,R^). The range of (g) 2 (-, 0 ) 
is equal to Qoo(Ao), so the minimal value, (g) 2 (xi( 0 ), 0 ) = inf( 5 oo(Ao) is situated 
strictly below infAeJ\{Ao} inf Qoo(A). 

In this paper, we shall work under the assumption that the subprincipal symbol of 
the unperturbed operator vanishes. 


Pi,o{x,^) = 0 . 


(2.28) 


The following is the main result of this work. 

Theorem 2.1 We adopt the assumptions above, in particular, 02.18|) . 02.211) . 02.27p . 
and 02.281) . Let us put x 1 (^ 2 ) = xi{f {^ 2 ), ^ 2 ) ■ Let 6 G (1/18,1/9) be fixed and assume 
that 


hWl-S) 


(2.29) 


Set 



11 

















Then for every Co > 0, we have the following description of the eigenvalues of P;. in 
the region 


\ h Tin z hi 

IzeC- |Re^|<-^, -<infQoo(Ao) + C'o^ (2.30) 

[ Oov^ ^ J 

valid for all h > 0 small enough: the eigenvalues are simple and given by 

\,k =p{f{h{j - O 2 )), h{j - 62 )) + ie{q) 2 {xi{h{j - 62 )), f{h{j - 62 )), h{j - 62 )) 

+ y/eh{X^ji. + + ...), 

(2.31) 

with j e Z, h{j — 62 ) = 0{h/y/e), N 3 k < 0(1), where Xf^ = Xl{h{j — 6 ^ 2 ), i/e) is 
a smooth function of ^2 = h{j — ^ 2 ) £ neigh(0, R) and \/e G neigh(0, R+), and 

A 2 ( 6 . 0 ) = e‘’'''‘(4p(/(6).&))’(C{9)2(^i(&)./(&).&))^ (*^+ 1 ) ■ (2'32) 

Here we have written 6*2 = fco(a 2)/4 + S 2 / 2 TTh, where ^ 0 ( 02 ) o,nd S 2 are the Maslov 
index and the classical action, respectively, of the fundamental cycle in Aq, given by 
a closed Hp-trajectory of minimal period. 

Remark. Choosing 5 G (1/18,1/9) in Theorem 2.1 to be close to 1/9, we see from 
fl2.29p that the description of the eigenvalues in Theorem 12.11 in the region fl2.30p is 
valid in the range 

when r; > 0 is small. In particular, we are able to reach some cases when e S> h, 
and here the analyticity assumptions seem essential. 

Remark. The result of Theorem 2.1 admits a natural extension to the case when 
Re^; G neigh(0,R) varies in a sufficiently small but hxed neighborhood of 0 G R. 
Indeed, let us recall the family of rational Lagrangian tori Ke C p~^{E), E G 
neigh(0, R), introduced in fl2.25p . A natural analog of the assumption fl2.2ip is then 
valid for infQoo(A£;), relative to the Lagrangian foliation in p~^{E), provided that 
\E\ is small enough. It follows therefore from Theorem 12.11 that the description 
fl2.3ip of the spectrum of remains valid when 

|Rez - E\ < < inf Qoo(Ae) + 
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uniformly in G neigh(0, R). We conclude therefore that the result of Theorem 12.II 
extends to the spectral region 


2 ; G C; |Re 2 ;| < 


1 Imz 


< inf Qoo(ARe 2 ) + o 




for C > 1 large enough. 


The plan of the paper is as follows. Section 3 is devoted to a general outline of the 
proof. In Section 4 we construct a global compactly supported weight function G, 
such that the leading symbol of P^, acting on the weighted space associated to G, 
becomes ~ p+ie {q — HpG), with the imaginary part avoiding the value £ inf Qoo(Ao) 
on p“^(0), away from the rational torus Aq. This effectively microlocalizes the 
spectral problem for to a small neighborhood of Aq. The quantum normal form 
construction for P^ in the rational region is carried out in Section 5, using the 
techniques of secular perturbation theory, thereby reducing the analysis to the study 
of a one-parameter family of non-selfadjoint operators in dimension one, having 
double characteristics with elliptic quadratic approximations. In Section 6 we recall 
the computation of low-lying eigenvalues for such operators, following |3] and |B], 
and extend the results there to the parameter-dependent case. The hnal step in the 
proof of Theorem 2.1 is taken in Section 7, where we carry out a pseudospectral 
analysis for the family of the one-dimensional operators in question, controlling the 
resolvent norms and obtaining the spectral localization. It then becomes possible 
to complete the proof by solving a suitable globally well-posed Grushin problem 
for Pe in a weighted space, using the ideas and techniques of 0.1121. In Section 
8 we present the results of numerical computations illustrating Theorem 2.1. The 
Appendix establishes some subelliptic resolvent bounds for non-selfadjoint operators 
of Schrodinger type, playing a principal role in the pseudospectral analysis of Section 
7 in the main text. These bounds seem to be of some independent interest, and their 
proofs are very much based on the techniques developed in a, 0. 


3 Outline of the proof 

In this section we shall give a general outline of the proof of Theorem 2.1. Some of 
the techniques come from the previous works [12], [10], and the presentation below 
will naturally focus on the new difficulties of pseudospectral nature, encountered in 
the analysis in the rational region. We shall then also describe heuristically some of 
the essential ideas employed in overcoming those difficulties, referring to Section 7 
and to the Appendix for a detailed rigorous discussion. 
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The principal symbol of the operator in fl2.ip . fl2.6p is of the form 


Pe=p + ieq + 0{e‘^), (3.1) 

in a neighborhood of p~^(0) r\T*M, and thanks to the ellipticity assumptions fl2.5p . 
fl2.9p . we observe that it suffices to make a microlocal study in the region where 
p is small. Recalling the assumption fl2.2ip and replacing q hj q — inf(5oo(Ao), 
in the following discussion we shall assume, for notational simplicity only, that 
inf Qoo(Ao) = 0. The hrst step in the argument is a construction of a global weight 
function G G C^{T*M) such that away from a small neighborhood of the rational 
Lagrangian torus Aq in p“^(0) fl T*M, we have 

(3.2) 


Away from Aq, the weight G satishes 

HpG = q — {q)j’i 

where {q)rp has been introduced in (12.141) . and when constructing G in a neighbor¬ 
hood of Aq, we introduce action-angle coordinates G T*T^, so that Aq = {.^ = 
0} C T*T2, and 

Pe{x, 0 = p(0 + 0 + C’(£^), (3.3) 

where the frequencies c}^^p(0) and d^^plO) are commensurable. After a linear change 
of variables, we get c?.c^p(0) = 0, and the isoenergetic condition (I2.18p shows that 
^liP(O) 7^ 0- following discussion, in order to hx the ideas, we shall consider 

the model case p(^) = ^2 + ^ 1 , which suffices to illustrate the difficulties. The weight 
function G near ^ = 0 satishes the cohomological equation 

HpG = q-q, (3.4) 

modulo C>(.^°°), where q = q{xi,^) is independent of X 2 and is such that 

g(xi,0) = — / q{x,0)dx2 (3.5) 

27r Jq 

is the average of q{x, 0) in the a; 2 -direction. From fl2.27l) we then know that q{xi, 0) > 
0 and that xi i-G- q{xi, 0) has a unique minimum which is nondegenerate. The partial 
Birkhoff normal form construction, utilized in solving fl3.4p may be continued, hrst 
at the principal symbol level, and then on the level of operators, leading to the 
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conclusion that microlocally in the rational region, when acting on an exponentially 
weighted space, the operator is unitarily equivalent to an operator of the form 

P{xi, hD^, E] h) + R{x, hD^, E] h) : L‘^{T^). (3.6) 

We refer to Proposition 7.1 in Section 7 for the precise statement. Here the full 

symbol of P{xi, hD^^^E] h) is independent of X 2 and is given by 

P(a;i, e; h) = p{^) + iEq{xi,^) + 0 {e‘^ + h‘^). (3.7) 

The contribution P(x, £; h) = 0{{e, h)°°) in fl3.6p is a remainder, which becomes 

0{h°°) when restricting the attention to the region ^ = 0{e^), for a suitable small 
hxed 5 > 0 — as we shall see, understanding this region suffices for the description 
of the eigenvalues in Theorem 2.1. In particular, since ^ becomes small, in the 
following heuristic discussion, we shall make a simplihcation and assume that q in 
fl3.7l) is independent of ^ altogether, depending on xi only. Let us also suppress the 
error term 0{e‘^ + h^) in fl3.7l) . for simplicity. When considering it in Section 7, it 
will be treated entirely as a perturbation. 

Taking a Fourier series decomposition in X 2 , we may view the operator P in fl3.6p as 
a one-parameter family of operators P{xi, hDxj^,^ 2 ,£] h) = P{^ 2 ), acting on L^(T), 
such that 

P{^2) = '^2 + -he, ^2 = hj, j E Z, (3.8) 

where 

+ iEq{xi), q>0, (3.9) 

is a one-dimensional non-selfadjoint Schrodinger operator with Eq as a potential. 
We are interested in the spectrum of the family 03.81) in the region where Rez is 
small and |Imz| < 0{h^/£), and the hrst observation is that the eigenvalues of the 
operator 

1 \ 2 _ ~ h 

-Le=[hD^^j +iq{xi), h =—, 

can be determined asymptotically in any disc |tc| < Ch, by means of the harmonic 
approximation, provided that h 1. See a, 0, and the discussion in Section 6 
below. The eigenvalues of in this region are of the form 

CXD 

hfc(h) ~ h e N, (3.10) 

j=0 
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where 


(3.11) 


I'i.o = h+ij, 

are the eigenvalues of the globally elliptic quadratic operator 

acting on L^(R). Here G T is the unique point such that = 0. The 

corresponding eigenvalues of P{^ 2 ) in (I3.8p are given by ^2 + £fJ>k{,h), and from [1], [6] 
we also know that 

II (^’(&) - WciLV) < O , (3.12) 

provided that |^ — ^ 2 ! < Chy/e and that (z — ^ 2 )/hy/e avoids the quadratic eigen¬ 
values fikfl in flS.lip . 

Now fl3.8l) is only an approximate direct sum decomposition, and in order to be able 
to absorb the error terms there, when constructing the resolvent of globally, it 
is of crucial importance to control the resolvent norms of also near the interval 
[Chy/e, 1/Ci(l)). To get such a control, since the spectral parameter remains close 
to the boundary of the range of the symbol of L^, we apply the method of ’’bounded 
exponential weights”, which in effect consists of replacing by a new operator 
for which the inhmum of the imaginary part is increased in the non-elliptic region 
for Re — oj). This method has been carried out in closely related situations 
in 13. a. 0 . and we apply some of those works in the actual proof in the Appendix. 
Here we shall merely recall the essential ideas. See also Ha, 122]. 

Let G(a;i,(^i) G be real-valued and odd in ,^ 1 . Let us consider formally the 
conjugated operator 

— ^-eG{xi,hD^^)/h o O Q^G{xi,hD^j^)lh^ 

acting on or equivalently, the operator acting on the weighted Hilbert space 
^eG{xx,hDx^)/hj^2^ y^^ Want this space to be equal to with its norm 

\\e-^G{xy,hD^^)/h^^y 

uniformly equivalent to the standard L^-norm. This is the case if the weight function 
G satishes suitable symbol estimates and has the fundamental property 

= Oil), (3.13) 
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uniformly with respect to the various parameters involved. 

We view e-^G{xT_,hDx^)/h ^ Pourier integral operator with the associated canonical 
transformation exp {ieHc), approximately equal to (xi, .^i) i—)■ {xi,^i)+ieHc{xi,^i), 
since £ will be small. Here Hq = ■ dx-^ — is the Hamilton vector held of 

G. By Egorov’s theorem we expect Lg to be an h-pseudodifferential operator with 
the symbol 

^ Le (exp {isHcixi^ii))) ^ L^{{xi,ii) + ieHG{L^){xi,ii)) 

Here Le(xi,.^i) = + ieq{xi) is the symbol of in fl3.9p . With ^{xl,^l) = we 

get 

Le{xi,ii) ^ + ie{q - Hi{G)){xi,ii) =: + ieg(xi, ^i). 

When considering — u for Rea; > h^/e, the most critical region is the one where 
^1^ ~ Rea; and it is here that we want to increase inf^.^ g as much as possible. 
Naturally, that will not be enough for the complete analysis, but in the following 
heuristic discussion, we shall restrict the attention to the region where = Rea;. 
Here, we get 


Qixi,^i) =q{xi) - 2^i4iG(xi,^i) 

=g(xi) — 2'/R^dxiG{xi, (Reo;)^/^), 

where we recall that G is odd in ^i, so that g is even in the same variable. Then 

zvReo; 

omitting = VRea; in the argument of G. We want 


inf g — inf g 

Xl Xl 


7^ 


(3.14) 


for a suitable small parameter 7, that we wish to have as large as possible, and to 
achieve this, we clearly have to modify g in a 7-neighborhood of x™“. Since we also 
wish IG*! to be as small as possible, we require 

supp G C [x“‘° - 7 , x““ + 7], 


and it is not hard to see that we can End such a G with 

„2 \ / »,3 


dx,G = O 


T 




u 


G = 0 


( 


\ 


u 
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The condition fl3.13p is fulfilled, provided that 




h\/Rea; 

Let C 3> 1 and let us choose 


= 0 ( 1 ) ^7 = 0 ( 1 ) 


hs (Rea;)6 


1 

53 


7 = — mm 
O 


1, 


hs (ReO') 6 


1 

53 


(3.15) 


(3.16) 


It follows from the heuristic discussion above that in the region where h^/e < Re a; < 
1/0(1) we obtain the spectral gain 


g7 

0 ( 1 ) 


mm 


in /i3 (Rea;)3£3 j > hy/e, 


in the sense that the resolvent (L^ — a;) ^ is well dehned in the region 

h^/e < Rea; < ^ , Imo; < -^min 73(Rea;)3£3,e'j 

C7 1 j (_/ V / 

and that in such a region, we have 

0 ( 1 ) 


(Lg — a;) ||/:(l2,l2) < 


min (hses (Rea;)3, e 


(3.17) 


(3.18) 


The resolvent estimates such as fl3.18p are established in the Appendix, using the 
machinery of bounded exponential weights and relying on the techniques of H. i- 
see Propositions A.2 and A.4 there, in particular. With the bounds fl3.18p available, 
we get the corresponding pseudospectral control over the family P{xi, hDx^,^ 2 , h) 
in (13.8p . in the region where |Re2; — ^2! > Chy/e, Im^; < 0{hy/e), and this allows 
us, eventually, to construct the resolvent of globally in this region. We therefore 
obtain some crucial spectral localization, making it possible to carry out the spec¬ 
tral analysis of working with one quantum number ,^2 = hj at a time, roughly 
speaking. A globally well-posed Grushin problem for P^ is hnally built from the 
corresponding one-dimensional Grushin problems for the operator in (13.911 . and 
solving it along the same lines as in 0.113. H. we complete the proof of Theorem 
2 . 1 . 

Remark. Our heuristic arguments seem to indicate that the optimal range for the 
perturbation parameter £ could be 

< 0^, (3.19) 
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as we need h = hjy/e 1 and e <^h. Due to many technicalities, we get a smaller 
range of values around £ ~ h, and leave extension to the range fl3.19p as an open 
problem for future works. 


4 Secular reduction and the global weight 

The purpose of this section is to construct a globally dehned compactly supported 
weight function, which will allow us to microlocalize the spectral problem for 
to a small neighborhood of the rational torus Aq. In doing so, we shall proceed 
similarly to |T2], with the essential difference that when working near the torus, 
the basic cohomological equation will have quite different properties, compared to 
the Diophantine analysis of [12], and will be treated using the secular perturbation 

theory, see ra. uni- 

Let us keep all the assumptions of Section 2 and consider the operator Pg with the 
leading symbol in fl2.10l) . in a neighborhood of p~^(0) fl T*M. Let 

/s:o : neigh(Ao, T*M) ^ neigh(e = 0, T*T^), (4.1) 

be a real analytic canonical transformation, given by the action-angle variables, such 
that the properties fl2.23l) hold. By Taylor expansion and fl2.23p . we have 

p(0 =P(/(6),6)+^(0 (Cl-/(6))% ^(0)>0, (4.2) 

where / is the analytic function introduced in fl2.24p . 

Implementing kq in fl4.1l) by means of a microlocally unitary multi-valued h-Fourier 
integral operator with a real phase, as explained in Theorem 2.4 in [7] and conjugat¬ 
ing Pg by this operator, we obtain a new h-pseudodifferential operator, still denoted 
by Pe, dehned microlocally near .^ = 0 in T*T^. The full symbol of P^ is holomorphic 
in a hxed complex neighborhood of .^ = 0, and the leading symbol is given by 

Ps{x, 0 = PiO + c) + 0{e^), (4.3) 

with p{^) of the form (14.21) . The function q in (14. 3 p is real on the real domain. On 
the operator level, Pg acts on the space of microlocally dehned Floquet periodic 
functions on T^, Lg(T^) C LjQp(R^), elements u of which satisfy 

O / 

u(x — u) = 6 = — - H —z/ e 27rZ^. (4.4) 

2'Kh 4 
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Here S = {Si, S 2 ) is given by the classical actions, 




Vdy, j = l,2, 


with aj forming a system of fundamental cycles in Aq, such that 

KQ{aj) = Pj, j = 1, 2, = {x e T^; X3_j = 0}. 

The tuple ko = (fco(«i), ^o(« 2 )) £ stands for the Maslov indices of the cycles aj, 

J = l,2. 

Remark. Using (14.2p . we see, using the implicit function theorem, that the energy 
surface p{^) = E, for E G neigh(0, R), is given by 


i2 + i{iuE) = 0, 
where £ is analytic with l{S,i, 0) ~ U^(0, 0) < 0. 


(4.5) 


Working near the zero section ^ = 0 in T*T^ and following the method of normal 
forms na, m, we shall now discuss the cohomological equation 

HpG = q-q, (4.6) 

where we want the remainder q to be simpler than q. Here we have 

Hp • dx, 

and thus, fl4.6p can be written more explicitly as follows, 

+ d^,p{^)dx,G = q-q. 

To simplify, we divide this equation by Writing u = G, v = {d^2P)~^Q^ ^ — 

we get 

{8x2 + a{^)dxi)u = V - V, (4-7) 

where a{^) = 9giP(0/^?2P(0- To simplify further, we replace the variables ^ by 


V= (Vl,V2) = (?1 - /(?2),&), 


(4.8) 


and write, abusing the notation slightly, u = u{x,rj), v = v{x,p), v = v{x,rj). It 
follows from (14■2p that the Taylor expansion of a has the form. 


a{7]) = ai{p 2 )yi + 02(^2)^! + . ■ ■, ai( 0 ) ^ 0 , 


(4.9) 
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and let us Taylor expand u, v and v similarly, 

OO 

= ^Uk{x,ri2)T]^, 
k=0 

OO 

= ^Vkix,T]2)vt (4.10) 

k=0 

OO 

= '^%{x,ri2)Vi- 

k=0 

Inserting these equations into (14. Th and identifying the powers of rji, we get 


d^^Uo = Vo-Vo, (4.11) 

d^^ui + aid^^uo = vi- vi, (4.12) 

9x2^2 + aidxiUi + a2dxiUo = V2 - T2, (4.13) 

and so on. The general equation is of the form 

dx^Uk + aidx^Uk-i + a2dx^Uk-2 + • • • + akdx^Uo = Vk - Vk- (4.14) 


u{x, rj) 
v{x, rj) 
v{x, Tj) 


The parameter 772 plays no essential role here and we sometimes suppress it from 
the notation. For a function u on the torus T^, we introduce its averages in Xk, 
k = 1,2, and its total average by 

1 

{u)k{x3-k) = — u{Xi,X2)dXk, 

271 Jq 

{{u)) = {{u)i)2 = 777 -^ / / u{xi, X2)dXidX2. 

[277) J 

Proposition 4.1 Let vq, vi,... G C'°°(T^) be smooth functions on T^. A neces¬ 
sary and sufficient condition on the smooth functions vq, G C°°(T^) for the 

existence of uq, mi, ... G C'°°(T^) solving fl4.11l) and fl4.14p for k > 1, is that 


{vo)2 = {vo)2, 
ii^k)) = {{vk)), k>l. 


(4.15) 


Proof: The necessity of fl4.15p follows from taking the a:2-mean of fl4.11l) and the 
total mean of fl4.14p . 
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Assume that the hrst equation in fl4.15p holds so that (vg — Vg )2 = 0. Then fl4.1ip 
has a solution ug = Uq ^ C°°(T^), given by 

rx2 

Ug{x)= {vg -Vg){Xi, t)dt. (4.16) 

Jo 

The general solution of (14.lip is of the form Uq{x) + /o(a:i), where fg{xi) is any 
smooth periodic function. 

We next consider fl4.12l) (i.e. fl4.14p with k = 1), which we write as 

dx2Ui = vi-vi- - aidxjg{xi). (4.17) 

Here the total average of vi —vi — aiOx^Ug vanishes, 

{{vi -vi- aidxiUg)2)i = 0 , 


and hence we can hnd a periodic smooth function fg{xi), unique up to a constant, 
such that 


(ui -vi- aidxiUg )2 - aidxjoixi) = 0. 

Equivalently, 

{vi -vi- aidx^Ug - aidxjg{xi ))2 = 0, 

and we can therefore hnd a solution G (^““(T^) to fl4.17p . and hence to fl4.12l) . 


Assume by induction that we have found ug, ui,... Uk-i, solving fl4.1ip and fl4.14|) . 
with k there replaced by j = 1, 2,..., /c — 1. We notice that the general solution 
of fl4.14p with k replaced by /c — 1, is of the form Uk-i = + fk-i{xi) for any 

smooth periodic function fk-i- We rewrite fl4.14p as 


dx2^k 'UJk Ojldxifk—l{pJl)i 


(4.18) 


where, 

Wk Vk Vk (^ldxi^k—1 2 

and we notice that {{wk) 2 )i = {{wk)) = 0. Choose fk-i such that {wk )2 = 
aidxj^fk-i{xi), or equivalently, so that {wk — aidxj^fk-i )2 = 0. Then there is a 
smooth periodic solution Uk = to 04.181) and hence to 04.14p . □ 

An application of Proposition 14.11 allows us to conclude that for any hxed iV G N, 
there exists an analytic function Go, dehned in a hxed neighborhood of ^ = 0, such 
that 

HpGo = q-q + 0((5i - /(&))~). (4.19) 
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for any analytic periodic fnnction q which satishes, 

( 9 (-, 0)2 = (g(-, 0 ) 2 , when = /(5) (4.20) 

and 

((?(•, 0)) = ((?(•> 0)), e neigh ( 0 , R^). (4.21) 

The following choice is convenient and will be made in what follows: let y : R ^ R 
be real analytic snch that x(0) = 0. We can then take 

fe, 0 = (1 - x (6 - /( 6 )))(?) 2 (xi, 0 + - /( 6 ))((g(-, 0)), (4-22) 

which is independent of X 2 - 

Let ns now restrict the attention to the energy snrface p“^(0). According to (14.51) . 
we have p(^) = 0 <;=> ^2 + -^(^i, 0) = 0, £(^ 1 , 0) ~ Q. It follows from (14.221) that when 
p (0 = 0 ) we may write 

inf g(xi, 0 = (1 - ^( 6 )) + ^^( 6 )^( 6 ), 


where 


fc(^i) = inf(g)2(a:i, 6, -^(6, 0)) = (g)2(2:i(6, -^(6, 0)), ^ 1 , -£(6, 0)), 

a:i 

9 i.i 1 ) — ((o'))('Cij 0 )), and -0 is an analytic fnnction snch that- 0 ( 0 ) = 0 ,'^'( 0 ) = 
x'(0). We compute next 

% inf g(a;i, 0 = ^'( 6 ) + '^'( 6 ) (^( 6 ) - + ' 0 ( 6 ) ( 6 ( 6 ) - ^'( 6 )) , 

XI 

dl inf 6 x 1,0 = fc"( 6 ) + 0 "( 6 ) ( 66 ) - 66 )) + 20 '( 6 ) ( 6 ( 6 ) - k'{i,)) 

^ tCi 

+ 0(6)(6'(6)-fc"(6)). 


Using that 


fc(0) = inf (g) 2 (a;i, 0 ) = infQoo(Ao) < (g)(Ao) = ((g))(0) = ^(0), 

xieT 

we see that the derivatives x'(0) and x"(0) of the analytic function y in (I4.22p can 
be chosen so that when p (0 = 0 , we have 

inf 6 x 1,0 > inf (g) 2 (xi, 0) + CO, (4.23) 

XI&T xieT 
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where the constant C > 0 is large. In other words, for A G neigh(Ao, J), we get 


inf q > inf Qoo(Ao) + 
A 


1 


dist(A, Ao)^. 


(4.24) 


Remark. In the preceding discnssion, we do not have to restrict onrselves to the 
energy snrface Indeed, introdncing the variables 

h = (hi,h2) = (6 - /(6),6), 

as in fl4.8p . and repeating the compntations above we get 

inf g(xi,0 > inf (9)2(^n/(6),6) + <^(6 -/(6))^ 

when p{^) = E, for E G neigh(0,R). Introdncing the rational Lagrangian tori 
Ag C dehned in fl2.25p . we therefore obtain on p~^{E), 

inf q{xi, C) > inf Qoo(Ae) + C'idist(A, KeY- (4-25) 

xiST 


We shall now constrnct a snitable global weight fnnction. In doing so, let Gt be an 
analytic fnnction dehned in a neighborhood of p~^(0) fl T*M, snch that 

HpCr = q — (?)r, (4.26) 

where 

1 

{q)^ = 7 f g o exp {tHp) dt, T > 0 
I J-T/2 

has been introdnced in fl2.14l) . An application of Lemma 2.4 of [12] together with 
the assnmption fl2.21l) allows us to conclude that outside an arbitrarily small neigh¬ 
borhood of Aq in p“^(0) n T*M, we have 

inf (g - L/pGr) > inf Qoo(Ao) + 7 ^, (4.27) 

<-'0 

provided that T is taken large enough. Here Gq > 0 is independent of the neighbor¬ 
hood taken. In these considerations, we are allowed to vary the real energy a little, 
and we conclude that for any hxed neighborhood W of 

IJ As, 0<Eo<l, (4.28) 

\E\<Eo 
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in p ^([—i?0) -^o]) there exists T large enough such that 


inf 

p-H[-Eo,Eo])\W 


(g — HpGx) > 


inf inf(5oo(A£;) + 

\E\<Eo Co 


(4.29) 


Here Co > 0 is independent of the neighborhood chosen. 

The global weight function will be obtained by gluing together the functions Gt ■= 
Gt o Kq ^ and Go, both viewed as analytic functions dehned in a neighborhood of 
the zero section = 0 in T*T^. Let -0 = G G°°(R^; [0,1]) be constant on each 
invariant torus ^ = Const, and assume that -0 = 1 near the rational region fl4.28|) . 
and with support in a small neighborhood of that set. Let us set. 


G = {1 - ij) Gt +'ipGo. 


(4.30) 


It follows that 

q-HpG = 'iPiq- HpGo) + (1 - ^)(g)r. (4.31) 

In a neighborhood of the rational region fl4.28p . we have 

q-H,G^q + 0 (Ki - /(&))") . 

with q given in fl4.22p . while further away from this set, we have q — HpG = {q)rp. 
In order to understand the behavior of {q)rp near ^ = 0 for T large, we write 

1 /■'^A 

{q)T{x,0 = 7f (l{x+ tp'{0,0dt, 

J- J-T/2 

and expanding g(-,0 ^ Fourier series, we obtain that 

(9>t(^.0= E e'^-^qlhOKiTk-p'iO)- (4.32) 

fc=(A:i,fc2)eZ2 


Here K is the Fourier transform of the characteristic function K of the interval 
[—1/2,1/2]. Let us decompose, 

(9>t(^.«) = E ■ k)q(k,i)e“-'‘ + E K(Tp'(^) ■ k)q(k,C)e“* = I + II. 

^27^0 k2=0 

(4.33) 

with the natural dehnitions of I and H. When estimating I, we use 04.21) and notice 
that when /c2 7^ 0, we have 

|p'(0 ■ > \p'iM - 0(1) I5i - /(6)l |4il > 1 - C |5i - /(&)! |4|. C > 0. 
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Here for notational simplicity we assume that the derivative of the function ^2 ^ 
p(/(^ 2 ),^ 2 ) is > 1 near 0. It follows that 

|P'(0 

provided that 2C |.^i — /(^ 2 )| 1^1 < 1- Let now 0 < y G C'^((—1,1)) be such that 
X = 1 on [—1/2,1/2] and write, 

I = ^x(2C|6-/(6)ll*:|)AWK)''=)».0e“-'' (4.34) 

fc2^0 

+ - x(2C-|6 - /(fe)| |fc|))/?(rp'(f) ■ k)q(k,Oe“* 

k2^0 

= 5: X(2C'I6 -/(6)l l*:|)o 

^2 7^0 

+ 5^(1 - x(2C|5i - /(&)! |*;|))A-(rp'(0 . fc)g(fc,Oe“-*. 

^ 2^0 

It is now easy to see, using the smoothness of q, that 

I = o(l + l 6 -/(fe)r). r>l. ( 4 . 35 ) 

When considering the contribution coming from II, we notice that 

n=(9)2(a:i,0+ E (A'(rpjyi)-l)e“''='?(i,0. (4.36) 

k2=0,ki^0 

Here | ~ |^i — /(^ 2 )|, in view of fl4.2p . and we conclude that in the rational region 
where = /(6), 6 e neigh(0,R), we get 

. 

Away from the rational region = f{^ 2 ), we see directly from fl4.33p that H con¬ 
verges to the torus average {{q)){0^ as T —)■ cx). 

Combining the equations and estimates (I4.19p . fl4.24p . fl4.26p . fl4.29p . and fl4.30p . we 
may summarize the discussion above in the following proposition. 
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Proposition 4.2 Let us make the assumption fl2.2ip . LetGo be an analytic solution 
near ^ = 0 of the equation fl4.6p . with q being of the form 04.221) . modulo — 

f{^ 2 ))^), for some N fixed large enough. There exists a real-valued function G G 
G^{T*M) such that G = Gqo kq in a neighborhood of Aq, and such that away from 
a small neighborhood of Aq in the region p~^{[—Eo,Eo\), 0 < Eq 1, we have 

g - i/pG > inf Qoo(Ao) + 7 ^, Go > 0. (4.37) 

When A C p~^(0), A G neigh(Ao, J), we have furthermore 

inf (g - HpG) > inf Qoo(Ao) + ^dist(A, Aq)^. 

A C/ 


Associated to the weight function G dehned in Proposition 4.2, we shall now in¬ 
troduce a suitable small but globally dehned deformation of thej^al phase space 
T*M into the complex domain. When doing so, let G G G^{T*M) be an almost 
holomorphic extension of G, and let us set 

A,g = exp C T*M. (4.38) 

Here a is the complex symplectic (2,0)-form on T*M and is the Hamilton 

vector held of Re G computed with respect to the real symplectic form Im a on T*M. 
It follows that the manifold A^g is I-Lagrangian and being a small deformation of 
T*M, it is also R-symplectic, i.e. an IR-manifold. From [19] and [24], we recall the 
general relation 

valid to inhnite order along the real domain T*M. Here ieH^ stands for the real 
vector held in T*M, naturally associated to the complex (1,0) vector held. 


ieH^ = is 


i=i 


(dG d 

\dfjdxj 


dG d\ 

dxj dfj j ■ 


It follows that in the region where G is analytic, including a sufficiently small but 
hxed neighborhood of Aq, we have 


A,G = exp{zeHG){T*M), 

where we write G also for the holomorphic extension and recall that exp {ieHG) is 
a holomorphic canonical transformation. 


27 









Associated to the IR-manifold A^g is the microlocally exponentially weighted Hil¬ 
bert space H{A^g)j dehned using the FBI-Bargmann approach, by modifying the 
exponential weight on the FBI transform side. We refer to IZD], HH for the detailed 
dehnition of the space //(Agc) in the case when M = R^, and to |26] and the 
Appendix of [7] for the case when M is compact. Following [19], [20], |2^, let us 
introduce a microlocally unitary h-Fourier integral operator 

Ug : L\M) ^ H{A,g), (4.39) 

dehned microlocally near p~^(0) r\T*M and associated to a suitable canonical trans¬ 
formation 

Kg : neigh(p"^(0), T*M) neigh(p"^(0), A^g), 
such that Kg = exp {isHg) near Aq. It follows that the operator 


P, : H{A,g) ^ H{A,g) (4.40) 


is microlocally near p ^(0) unitarily equivalent to the conjugated operator 


Ug^P^Ug : ^ L\ 


with the leading symbol 


PeWa ^P + i£{q- HpG) 0{e‘^). 


Letting 


Uo : L\M) ^ 


be the semiclassical microlocally unitary Fourier integral operator with a real phase 
associated to the canonical transformation kq in fl4.1l) and using the operator UqUg^ 
associated to the canonical transformation 


Ko o Kg^ : neigh(exp {ieHG){Ao), A^g) neigh(^ = 0,T*T^), 

we get that microlocally near the Lagrangian torus exp (i£iLG)(Ao) C A^g, the 
operator in fl4.40p is unitarily equivalent to an operator P^, acting on Lg(T^), dehned 
microlocally near ^ = 0 in T*T^, given by 

OO 

R, ~ (4.41) 

u=0 
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Here are holomorphic functions in a fixed complex neighborhood of ^ = 0, smooth 
in £ G neigh(0, R), and 

Po = P(0 + 0{e^)+eO{{^i - /(6))^), (4.42) 


with q{xi,^) independent of X 2 and of the form fl4.22p . Furthermore, the assumption 
fl2.28p implies that 

= 0{e). 


We may illustrate the microlocal unitary equivalence above by the following com¬ 
mutative diagram, 


Pe : H{A,g) 


H{Kg) 


UoU-. 


Pe : 


UoUr, 




(4.43) 


In what follows, we shall drop the tildes from the notation in fl4.4ip and write simply 
Pe and py, 1 / > 0. 


5 Quantum normal forms near rational tori 

In this section, we shall be concerned with a classical h-pseudodifferential operator 
Pe {x, hDx] h), dehned microlocally near .^ = 0 in given by the expansion 

fl4.4ip . with the leading symbol of the form fl4.42p . Our purpose here is to obtain 
a normal secular reduction of P^, also on the level of lower order symbols, and this 
will be accomplished in a way very similar to n, Hi. 

Let us hrst discuss the normal form construction at the level of principal symbols. 
In doing so, we let % := q in 04.421) . and write 

Po{x, e) = p(0 + ieqo{xi,^) + ie‘^qi{x, 0 + 0 + £(6 - /(6))^) • (5-1) 

Arguing as in Section 3, we can construct an analytic function Gi, dehned near 
.^ = 0, such that modulo O{{0 ~ we have 

HpGi = qi — qi-, 

where qi is any analytic function satisfying 04.201) . 04.2ip . with q replaced there by 
qi- It follows that 

Po(exp {ie^HG^){x,0) = p(0 + 0 + 0 + C»(£^ + £(^i -/(6))^). 

Continuing this procedure, we get the following result. 
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Proposition 5.1 Let po{x, e) = p{^) + ieqo{xi,^) + Oie"^) + eO{{^i - /(6))^) be 
analytic defined near ^ = 0, depending smoothly on e E neigh(0,R). Here N > 2 is 
arbitrarily large but fixed. Assume that 

pio = p(/(e 2 ),6) + gm > o, /(o) = o, 

where p{f{^ 2 ),^ 2 ) = C (^2 + « > 0. There exists a holomorphic canonical 

transformation of the form 

kW = exp o ... o exp (5.2) 

with Gj analytic near ^ = 0, 1 < j < iV — 1, such that modulo an error term of the 
form + £(^i — f{f 2 ))^), we have 

Po{nf\x, 0) = p(0 + ieQoixi, f) + ie^{qi{xi,f) + ... + 

is independent of X 2 . Here, as discussed before, go is any analytic function satisfying 
fl4.20p . fl4.2ip . and inductively (fk is any analytic function satisfying fl4.20p . fl4.2ip . 
with q there replaced by a certain function qu that depends on the previously chosen 
%■,■■■ Qk-i- 


As will be discussed in Section 7, the complex canonical transformation in 
fl5.2p can be quantized by means of an elliptic classical h-Fourier integral operator 
74 in the complex domain, depending smoothly on e: G neigh(0,R), introduced 
rigorously on the FBI transform side. In this section, we shall proceed formally, and 
an application of Egorov’s theorem allows us to conclude that the operator 

Pfix, hD,,] h) = Ufi^Pe{x, hD,c] h)Ue 

is an h-pseudodifferential operator, dehned microlocally near = 0, whose symbol 
has a complete asymptotic expansion 

Pe{.x, h) ~ Po + + ..., (5.3) 

with all pj = Pj{x, f, e) being smooth functions of £ G neigh(0, R) with values in the 
space of holomorphic functions in a hxed complex neighborhood of ^ = 0, such that 

Po{x, 4 e) = p{i) + ieqo{xi,i) + 0{e^) + + £(^i - /(6))^). (5.4) 

Here the Ci(£^)-term is independent of X 2 and has the properties described in Propo¬ 
sition EUl Furthermore, we still have pi{x,f,e) = 0{e). 
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We shall now simplify the lower order terms Pj, j > 1, in fl5.3p . To that end, let 
y4e(x, hD] h) be a classical analytic elliptic h-pseudodifferential operator of order 0, 
with symbol 

A^{x, h) ~ ao{x, e) + hai{x, ^,e) + (5.5) 

depending smoothly on e: G neigh (0, R). Then 

A;^U-^P,U,A, = A-^P,A, =: P,(x, hD^- h), 


where 

Pe{x, h) ~ po{x, e) + hpi{x, e) + h^p 2 {x, e) + h^%{x, e) + ..., (5.6) 

with ^ ^ 

Pi=Pi + -aQ^Hp^ao =pi + -Hp^bo, (5.7) 

if bo = Inoo, well-dehned up to a constant. Thus, looking for bo in terms of a formal 
power series in e and choosing the terms there suitably, we can arrange so that 

= Pi,o(xi, 0 + + • • • + + (^1 - /(6))^), 

where pi^o is any analytic function satisfying 04.201) . 04.2ip with q replaced by pi,£=o, 
and inductively, pi^k is any analytic function satisfying 04.201) . 04.2ip . with q replaced 
by a function depending on the previously chosen pi q, ...,pi^k-i- 

Iterating this procedure, by choosing also the lower order terms in the expansion of 
Ae, we get the following result, giving a quantum secular normal form construction. 


Proposition 5.2 Let 

Pe ^ Po +hpi + ... (x, 0 e neigh(^ = 0, T*T^), 

be such that po{x,^,e) has the properties stated in Proposition 15.11 Let he an 
elliptic classical analytic h-Fourier integral operator of order 0, associated to the 
canonical transformation in Proposition \5A[ Then we can construct an elliptic 
classical analytic h-pseudodifferential operator of order 0 with symbol as in 05.5p . 
such that 

A;^U-^P,{x, hD- h)U,A, = P,{x, hD^-h), (5.8) 

where P^ix, hD^, h) is of the form 

Pe{x, h) ~ po{x, f, e) + hpi{x, f, e) + h^p 2 {x, f,e) + .... (5.9) 
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Here the leading term po is as in fl5.4p 


e) = p(0 + ieqo{xi,i) + 0 {e^) + + e(^i - /(6))^), 

with the 0{e‘^)-term being independent of X 2 - For 1 < k < N, we have 

Pk{x, f, e) = Pk,o{xi,0 + 0 + • • • + + (^1 - /(6))^), (5.10) 

where pk/ is any function satisfying fl4.20p . fl4.21l) . with q replaced bypk/; a function 
that depends on p-j^j, for all {k,i) such that either k < k, or k = k and i < i. In 

particular, we can choose pk independent of X 2 modulo + {fi — f{f 2 ))^)- We 

also have 

= 0 {e). 


6 Harmonic approximation for non-selfadjoint o- 
perators 

In the previous section, we have seen how to eliminate the 0 : 2 -dependence in the 
complete symbol of our operator, by means of successive averaging procedures, when 
working in a small neighborhood of the rational torus Aq = {.^ = 0} C T*T^. 
Following Proposition 15.21 and neglecting the remainder terms there, we shall now 
consider an operator of the form 

Pe = Pe{xi,hDa;-,h), (6.1) 

dehned microlocally near .^ = 0 in and acting on with a complete 

symbol independent of X 2 - We assume that the leading symbol of P^ is of the form 

PoOui.e) = p(0 + + C>{e^), p(0 = p(/(6), ft) + g(i) (ft - /(ft))^, 

( 6 . 2 ) 

with q given in 04.221) . and let us recall the assumption 02.27p implying that the 
function T 3 Xi 1 —)■ q{xi, f{f 2 ),^ 2 ) has a unique minimum which is non-degenerate, 
for ^2 G neigh(0,R). When discussing the spectral analysis of P^, it is natural, 
in view of its independence of X 2 , to take a Fourier series expansion in X 2 , thereby 
reducing the problem, at least formally, to a direct sum of one-dimensional operators 

p,{x,,hD,„h{j-e 2 y,h), jez, = + 

4 2 nh 

considered for those j for which h{j — 62 ) G neigh(0, R). 
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In what follows, we shall write ^2 = h{j — 62 ) G neigh(0,R), and concentrate our 
attention on the one-dimensional operator 


Pe{xi,hD^^,^ 2 -,h), (6.3) 

acting on Lg^(T). Modifying the Floquet conditions on T, we may replace (I6.3p by 
the conjugated operator 

The full symbol of Pe{xi, + /(^ 2 ), ^ 2 ; is of the form 


^£(2:1,6 + /(6), 6; h) = Y ^ h^pj^eixi,^). 

j=o 


Here 


PoA^i^ 0 = p(6 + /(6), 6) + i£q{xi,^i + /(6), 6) + 0 {e^), 

and from Proposition 15.21 we recall that 

piA^uO = c’(£). 


(6.4) 

(6.5) 

( 6 . 6 ) 


We can then write pi^e = £qi,e, qi,e = 0{1). 


Let us set 


h 

(6.7) 

h=^, 

Ve 

and assume that 


h < 1. 

(6.8) 

We have 



Pe{Xu^l + f{^2),^2;h) 

00 

= P(6 + /(6), 6) + ieq{xi,^i + /(6), 6) + + heqi,e{xi, 0 + h^P. 


i=2 


p{^i + f{^2),^2)+£ I iq{xi,^i + f{^2),^2) + 0{e) + he^Pqi^e + 

i=2 


Here, according to fl6.2p . 


p(6 + /(6),6) = p(/(6),6) + ^7(6 + /(6),6)e?, 77(0) > 0. 
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It follows that on the operator level, we have 


Peixi, /fe) + hD.i, 6; h) = p(/(6), 6) + 17(/(6) + f 

+ ie (^q{xi, f{^ 2 ) + hD^-^,^ 2 ) + 0{e + he^P + . (6.9) 

We shall be interested in computing eigenvalues of the one-dimensional operator 
Pe{xi, /(^ 2 ) + hD^^,^ 2 ; h) in the region 

\Rez-p{f{^2),^2) \ < 0{€h), 


and directly from fl6.9p . using cut-offs of the form x(^i/\/£), we see that the cor¬ 
responding eigenfunctions are microlocally concentrated to the region where ^1 = 
0 {^/e), provided that the smallness condition fl6.8p is strengthened to the following 
one, 



?7 > 0. 


( 6 . 10 ) 


It will then be convenient to perform a rescaling of the cotangent variable, corre¬ 
sponding to a suitable change of the semiclassical parameter. Let us write 


hPxi x\ ; 


and if ^ 1 , denote the cotangent variables corresponding to and hD^^, re¬ 
spectively, we have ^ 

6 = 

It follows that 

^PeixiJ{^2) + hD^^,^2;h) ( 6 . 11 ) 

can be viewed as an h-pseudodifferential operator of the form 


-Pe{xiJ{^2) + hD^^,^2]h) = 
e 

+ ,^ 2 ) 


p(/(6),6) 

£ 

{hD^^ + iq{xi, f{^2) + x/eh/lxi, 6) + C>(£) 


+ hO{y/E) + 0{h^). (6.12) 


Ignoring the constant term p{f {^ 2 ), ^ 2 ) /£ in the right hand side, we recognize here es¬ 
sentially a one-dimensional Schrodinger operator with a purely imaginary potential, 
and to be precise, we can write 

--^£(2^1; /('C 2 ) + h) = pHIAJjAJ. _|_ A[xi, hDx^,^2, V^; h), 
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where A{xi,hDxj^,^2,'/£]h,) is a well-behaved h-pseudodifferential operator, de¬ 
pending smoothly on ^2 ^ neigh(0, R) and y/e > 0, with the leading symbol 

gif {^ 2 } + f (^ 2 ) + y/^^ 1 ,^ 2 ) + Oie), (6.13) 

and with a subprincipal symbol which is Oiy/e). Here we have dropped the tilde 
from the notation for the cotangent variable corresponding to and let us also 

recall that the operator fl6.12l) is to be considered microlocally in the region where 
= 0(1). The function g in fl6.13p satishes g > 0. 

Remark. If fl6.6p is no longer assumed, we can write, assuming that h/e <C 1, 

-^£(2:1, + /(C2), C2; h) = pi^i -\- 7(^2), ^2) + Cl + /(C2), C2) 

00 

-h 0(£^) -F hpi^eixi, C) + X] h^Vj,ei^l,0 

1=2 

= p(Ci + /(C2), 6 ) + ie ^g(xi,Cl + /(C2),C2) + o(^ + ^ + ^ 

and we can then view h/e as an additional small parameter. As will be seen in 
Section 6, some pseudospectral considerations will force us to assume that e/h should 
not be too large, and for that reason, in this work we make the assumption fl2.28|) . 
leading to fl6.6p . 

The discussion pursued in this section so far indicates that the spectral analysis of 
the original operator should reduce to that for a family of h-pseudodifferential 
operators on T, with leading symbols of the form (16.131) . Letting e = 0 in (16.131) for 
a while and suppressing the parameter C 2 altogether, we shall now pause to make a 
digression, in order to recall semiclassical asymptotics for the low lying eigenvalues of 
non-selfadjoint h-pseudodifferential operators with double characteristics. In doing 
so, we shall follow the analysis of [B], which in turn follows [1] closely. Let us also 
remark that in the present one-dimensional case, the quadratic approximations along 
the double characteristics are elliptic and consequently, our discussion is considerably 
simplihed, when compared with i, i. 

Let P„{x,hD,;h) : C“(T) ^ C“(T) be such that P„(x,hD,;h) € Op”(S({0^)), 
and assume that the semiclassical leading symbol of Pq is of the form 

Poix,0 = + '^yix), (6.14) 
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where V E C'°°(T;R). Assume also, for simplicity, that the subprincipal symbol of 
Pq vanishes. Let us assume that if a = mini/ then V~^(a) = {xq} with V"(xo) > 0. 
We are interested in the eigenvalues of Pq in an open disc {z E C-,\z — ia\ < Ch}, 
for some C > 0 hxed and all h > 0 small enough, and to that end we consider the 
operator 

P{x, hD^] h) = {1 — i){PQ{x, hD^i h) — ia), (6.15) 

whose leading symbol p{x,^) = (1 — i){po{x,^) — ia) is such that 

Rep{x, 0 = + ^( 2 ^) “ o > 0 

is elliptic for large ^ and vanishes precisely at the point (xq, 0) G T*T. In a neigh¬ 
borhood of (xo, 0) we have 

p(x-/xo,0 = ?(a^,0 + C’((x,0^), (a;,O^(0,0), (6.16) 

where g is a quadratic form, such that Reg > 0. When determining the eigenvalues 
of P{x, hDx] h) in an (9(h)-neighborhood of 0, naturally only the behavior of the 
operator in a small neighborhood of (xq, 0) matters, and by composing p with an 
inverse of the translation 


K : neigh((xo,0),T*T)neigh((0,0),T*R), k((xo, 0)) = (0, 0), (6.17) 

we obtain an h-pseudodifferential operator 

00 

P{.x,i]h) ^^ypj{x,i), Pi = 0, (6.18) 

3=0 

dehned microlocally near (0, 0) G T*R, such that the leading symbol Po = p satishes 

P{x, 0 = Qix, 0 + 0^), (6-19) 


where g is quadratic with 


Reg > 0. (6.20) 

We extend P(x, h) to be globally dehned on R^ as an element of the symbol class 
S'(l), such that 

Rep(x, 0 > 0, (Rep)"^ (0) = {(0, 0)}, (6.21) 


and such that 

Rep(x,0>^, |(2:,0I > C'> 0. (6.22) 

An application of Theorem 1.1 of [6] allows us to conclude that the following result 
holds, which we state directly for the operator Po{x, hD^] h). See also [5] for related 
results in the analytic case. 
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Theorem 6.1 Let the operator Pq{x, hD^; h) : C'°°(T) —)■ C°°(T) have the principal 
symbol of the form fl6.14l) . a vanishing subprincipal symbol, and let us assume that if 
a = minV^ then V~^{a) = {a:o} with b := V''{xq) > 0. Let C > 0. Then there exists 
/lo > 0 such that for all 0 < h < ho, the spectrum of the operator Po{x, hD^] h) in 
the open disc in the complex plane D{ia,Ch) is given by the simple eigenvalues of 
the form, 

Zk ^ ia + h (Afc^o + hX^ i + h^\k ,2 + •••)• (6.23) 

Here Ak,o the eigenvalues in D{f],C) of the elliptic guadratic operator 

= Dl + i^x"^, 

acting on L^(R), which are given by 

1/2 

Ko=[^] e^^/\2k + l), fcGN, k = 0{l). 

Remark. Theorem 16.11 continues to be valid when the operator Pq{x, hD^] h) acts 
on an L^-space of Floquet periodic functions on T and indeed, the eigenvalues 
described in this result do not depend on the Floquet conditions, modulo 0{h°°). 

Coming back to the operator in fl6.12l) . with the leading symbol (16.131) . we shall 
next have to extend the result of Theorem 16.11 to the parameter dependent case, 
and to this end it will be convenient to recall briefly the main steps in the proof 
of Theorem 16.11 Let P = P{x, hD^] h) be an h-pseudodifferential operator on R 
satisfying fl6.18|) - (16.221) . Following [6] , let us recall that the proof of Theorem 16.11 
proceeds by constructing a well-posed Grushin problem for the operator P, of the 
form 

{P — hz)u + R-U-= V, R+u = v, z G neigh(Ao, C), (6.24) 

in the space T^(R) x C. Here Aq is an eigenvalue of q^{x,Dx) such that |Ao| < C. 
The operators P_ : C —)■ and P+ : —)■ C are dehned as follows, 

R^u_ = u_e, Rj^u = {u, f)i 2 , (6.25) 

where e is an eigenfunction of q'^{x, hD^) corresponding to the eigenvalue hAo, and / 
is an eigenfunction of the adjoint operator q{x, hD^), corresponding to the eigenvalue 
hXo. 

The verihcation of the well-posedness of (I6.24|) consists of two steps, both carried 
out after a metaplectic FBI transform, 

T : L‘^{K) ^ H^,{C). (6.26) 
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Here 


H^^{C) = Hol(C) n L\C; e-^^°/’^L{dx)), 

and $0 is a suitable strictly subharmonic quadratic form. In the first step, we 
concentrate on the region |x| < h^, x G C, for some 1/3 < p < 1/2. Arguing as 
in [6], we obtain the following a priori estimate for the problem (I6.24p . based on the 
quadratic approximation of P near the origin — see formula (3.25) in [6], 

\\{h+ (^) “II + l“-l 

< C\\{h+\xf)-^^\o n|| + C|| {h+\xf)-^^\o {P-Q)u\\ 

+ |n+| + {h + \xfy/^lK (^) u\\. (6.27) 

Here u, v E i7<i,g(C), the norms are taken in the space L^{C; e~^^°^^L{dx)), and we 
have also written P for the conjugated operator TPT~^. The function xo ^ C'^(C) 
is fixed, with Xo = 1 near 0, and A' is a fixed compact neighborhood of supp(Vxo), 
0 ^ A". Furthermore, Q = Tq'^{x,hDx)T~^, and therefore, as explained in [6], we 
have 

II (A + (|;) (P -Q)u\\ = 0 II (6.28) 

Using fl6.27p and (16.281) . we obtain 

A||Xo (^) “It+ l“-t 

< ^11« II" + 0(A«''-‘)|| u II" + 0(h) |t,+ |" + 0(A)|| 1*- (|;) ti II". (6.29) 

Notice that the lower bound p > 1/3 implies that here h. 


In the second step of the proof, we consider the exterior region, |x| > and here 
we use that 


Rep 


X, 




dx 




U > 0. 


Exploiting the sharp Carding inequality in the form of a quantization-multiplication 
formula, as explained in [^, see also [29], we obtain the following exterior a priori 
estimate for the problem (16.241) . 





\u{x) 


g-2$o(UA < 0(1)11 n II II M 


+ 0(h°°)|M_| ||n||+ 0(h)||M||2 


(6.30) 
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Here x ^ [0,1]) vanishes near a; = 0 and x = 1 large x. Assnming that 

1/3 < p < 1/2, the bonnds (16.291) and (16.301) can be glned together, and we get the 
a priori estimate 

h\\u\ \ + \u-\ < C>(1) (II n II + h |n+|), (6.31) 

and the conseqnent well-posedness of the Grnshin problem (16.241) . Asymptotic ex¬ 
pansions for the eigenvalnes of P follow exactly as explained in 01,0. In the present 
one-dimensional sitnation, the eigenvalnes are simple and only integer powers of h 
occnr in the expansions fl6.23p . 


Tnrning the attention to the parameter-dependent case, let Pe = Pe{x, hD^; h) G 
OPh(*S'(l)), e > 0, be an h-psendodifferential operator depending smoothly on y/e, 
snch that Ps=o = P satisfies flb.lSp . (16.191) . fl6.20p . fl6.2ip . fl6.22p . In particniar, the 
leading symbol of P^ satishes 


Pe{x,0 =P{x,0 (6.32) 

in the sense of symbols in S'(l), and the snbprincipal symbol of P^ is 0{y/e). Assnme 
also that near (0,0), fl6.32p improves to 

Pe{x,0=p{x,0 + (^{y^\^\+^)^ (6.33) 

see also fl6.13p . We wonld like to conclnde that the Grnshin problem fl6.24p with P 
replaced by P^ remains well-posed, provided that £ > 0 is not too large, and to that 
end, we shall simply inspect the two steps above. 

In the region |x| < h^, we argne as above, with P replaced by P^, and nsing fl6.33p . 
together with the fact that snbprincipal symbol of Ps is 0 {y/e), we see that we get 
an additional term in the right hand side of (I6.27p of the form 

11 (A + (I)) (P, - P)'« 11 = O 11« 11. (6.34) 

Here we also assnme that we have chosen the FBI transform in (I6.26p so that (I6.33p 
holds on the transform side. As for the exterior region |a:| > h^, replacing P by Pg, 
we get an additional term in the right hand side of (16.301) . given by 

0{l)^/e\\u\\^. (6.35) 


It follows that to absorb the two extra terms (I6.34p . (I6.35p . we need to meet the 
following conditions. 


hV2 


« 
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and 


The first condition is satisfied provided that 

since p < 1/2, and the second one holds when 

£ < h^P. 

We conclude that the Grushin problem fl6.24p remains well-posed when P is replaced 
by Pe, provided that 

e < (6.36) 

since 1/3 < p < 1/2. Combining this observation with the standard perturbation 
theory for eigenvalues of multiplicity one [13], we obtain the following result. 

Proposition 6.2 Let Pe{x, hD^, h), e > 0, be a smooth function of s/e with values 
in Op)('(S'(l)), such that when e = 0 , we have the properties fl6.18p - fl6.22p . Assume 
that fl6.33p holds. Then for e < , p > t), the eigenvalues of P^{x, hD^] h) in the 

region {z G C; \z\ < Ch} are given by the simple eigenvalues of the form 

Zk ^ h (Afc^o(V^) + hXk^i{\/e) fcGN, k = 0(1), 

where \kj{\fe) are smooth functions of s/e > 0, j > 0, with Afc^o(O) being the 
eigenvalues of the guadratic operator q'^{x, D/j, described explicitly in Theorem 16.11 
When z E C is such that \z\ < Ch and dist( 2 ;. Spec (P^)) > h/0{l), we have 

{z - P/~^ = ^ L\ (6.37) 


In our considerations, see fl6.13l) . when applying Proposition 16.21 we should replace 
the semiclassical parameter h hj h = h/s/e, which in view of fl6.36p leads to the 
condition ^ 

e<^h^P, l/3<p<l/2, (6.38) 


so that 


£ <C h 2p+i. 


When p = 1/3, the power in the right hand side = 4/5, and it follows that we have 
the well-posedness of the Grushin problem provided that 

e < 0(ht+^), r/ > 0. (6.39) 
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Remark. In the proof of Proposition 6.2 above, the presence of the parameter ^/e 
was treated by a direct perturbation argument, leading to the upper bound fl6.36p . 
The purpose of this remark is to outline an alternative approach to the parameter- 
dependent case, leading to sharper bounds on e. While sharpening the result of 
Proposition 6.2 below would not lead to an improvement in Theorem 2.1, which 
is the main result of this work, we believe that the alternative approach sketched 
below may be of some independent interest. Since its precise realization is likely to 
demand a greater technical investment, the argument developed in this remark will 
be quite brief and we hope to be able to develop it further in a future work. 

Let h) be a real analytic function of e G neigh(0, R) with values in the space 

of bounded holomorphic functions in a tubular neighborhood of R^, such that as 
h 0+, 

Pe{p] h) ~ p,{p) + + ..., p = (x, 0- 

For £ = 0, let us assume that the leading symbol p := po is such that Rep > 0 is 
elliptic at inhnity, vanishing precisely at p = 0. Assume furthermore that we have, 

P{p) = (l{p) + 0{p^), p^O, 

where q is quadratic with Reg positive dehnite. In particular, p = 0 is a non¬ 
degenerate critical point for p and an application of the implicit function theorem 
shows that for e small, p^ has a non-degenerate critical point p{e) in the complex 
domain, depending analytically on e, with p{e) = 0{e). Passing to the FBI trans¬ 
form side by means of a metaplectic FBI transform T, as in fl6.26p . let us continue 
to write p{e) = (a;(£),^(£)) = 0(e) for the image of the critical point p{e) under the 
complex linear canonical transformation kt associated to T. 

We know from m that kt(R^) = A$q = {(x, {2/i)dx^o(x)); x G C}, where <Fo is 
the strictly subharmonic quadratic form introduced in fl6.26p . We shall now discuss 
the problem of constructing a weight function G C'°°(C) such that 

$, = $0 + 0(h), IV' (<!>, - <I>o)I < 1, (6.40) 

and with p{e) G = {(x, {2/i)dx^£{x)); x G C}. The function $£ is then strictly 
subharmonic and if we set iL#^(C) = IIol(C) fl L‘^{C;e~‘^'^P^L{dx)), then we have 
with uniformly equivalent norms. To get the complete asymptotic 
expansions of the eigenvalues of in D{p^{p^),Ch), as in Proposition 6.2, one 
should then work with the operator P^ acting on the space . We need 
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and let us notice that 


m - = o(£)' 

With dx^£(x(e)) already determined, we try 
<he(a;) = <ho(x) + 2Re - aa;<l'o(a;(£))) • (x - x(e))^ x 

= %(x) + h“ (4x) 

Here y G C“(C) is a standard cut-off near 0 and 


X — x(e) 

X — x(e) 


4(l/) = 2Re - a 3 ,<ho(a;(£))) • y 

is linear, such that 4 = as a linear form. Then 

($, - $o) = k>o, 


and in view of (I6.40p . we need eh° < 0(h), e/h" -C 1. We get the conditions 
£ < e 0{h°'), and it follows that the optimal choice of a is given by 

a = 1/2. This leads to the condition e -C 0{'/h). In our applications, we should 
replace £ by and hhj h = h/^/e, leading to the condition 

T h 
e <^h = 

so that we get 

0 < e < 0{h^l^), (6.41) 

which is sharper than fl6.39p . One conjectures therefore that the result of Proposition 
6.2 extends to this range of e and we hope to return to this observation in a future 
paper. 


We shall hnish this section by a formal application of Proposition 16.21 to the mi- 
crolocally dehned operator Ps{xi,hDxj^,^ 2 ',h) in fl6.3p . acting on Lg^(T): assume 
that £ > 0 is such that 

h = T] > 0, 

Ve 


and that fl6.39p holds. It follows that the eigenvalues of Ps{xi,hDxj^,^ 2 '-,h) in the 
region 

4-p(/(6),6) -*^(?)2(a:i(6),/(6),6)l < 0 {y/eh) 
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are given by 


Zk =p(/(6),6) +^^(?)2(a;i(6),/(6),6) 

+ y/eh(Xk^Q + Xk^ih + + •••)) N 3 fc < 0(1), 

where = Xkj{^ 2 , V^), j > 0, is a smooth function of ^2 £ neigh(0,R), 
with 

Am« 2.0) =e"'-‘(S|.p(/K2).?2))i(C{9)2(^i(&)./(&).&))4 (*^+ ^) ■ 

Here we recall from fl2.27p that a;i(.^ 2 ) € T is the unique point of minimum of the 
function xi H- (g) 2 (a;i,/(6), 6)- 

7 Pseudospectral bounds and the global Grushin 
problem 

The discussion pursued in the previous section shows that we are able to determine 
the low-lying eigenvalues of suitable localized one-dimensional operators 

Peixi, hD 

XI 7 ■^ 2 ; h) 

in fl6.3p . occurring in the normal form reduction, provided that the perturbative 
parameter 5 satishes 

<e< r]>0. (7.1) 

The purpose of this section is to complete the proof of Theorem 12.11 by construct¬ 
ing a global well-posed Grushin problem for — z, leading to the description 
of the eigenvalues in the region described in Theorem 12.11 In doing so, we shall 
have to strengthen the bounds in fl7.ip . as a consequence of some precise pseu¬ 
dospectral analysis for the family of the one-dimensional non-selfadjoint operators 
Pe{xi, hDxi,^ 2 ] h), with ^2 playing the role of parameters. 

Our hrst task is to give a global dehnition of the h-dependent weighted Hilbert 
space, where the Grushin problem will be studied. Similarly to |T2], the weighted 
space in question will be associated to a globally dehned IR-manifold A C T*M, 
which is C)(e:)-close to T*M and agrees with it outside a compact set. Specihcally, 
the manifold A will be obtained as an (!l(£^)-perturbation of the IR-manifold A^g, 
introduced in fl4.38p . where the perturbative modihcation will only take place in a 
sufficiently small but hxed neighborhood of the rational torus Aq. 


(6.42) 
^/e > 0 , 

(6.43) 
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Let us recall therefore that in Section 4, we have shown that microlocally near 
the Lagrangian torus exp (ieiLG)(Ao) C A^g, the operator in fl4.4Up is unitarily 
equivalent to an analytic h-pseudodifferential operator dehned microlocally near 
.^ = 0 in and acting on Lg(T^), such that the leading symbol of Pg is of the 

form 

Po{x, e) = p(0 + m{xi, 0 + 0{e^) + eC>((6 - /(6))^), (7.2) 

where p{^) is given in fl4.2p and > 2 is arbitrarily large but hxed. See also fl4.43p 
for an illustration of the unitary equivalence by means of a commutative diagram. 

Let 


4^) : neigh(^ = 0, T*T‘^) neigh(^ = 0, T^T^), j21x7? (7.3) 

be the holomorphic canonical transformation, introduced in Proposition 15.11 Con¬ 
sidering the IR-manifold Ki^\T*T‘^) C T*T^, dehned in a complex neighborhood of 
^ = 0, we conclude, arguing as in Section 5 in [12], that there exists a strictly 
plurisubharmonic function $£( 0 ;), dehned for x G |Imx| < 1/(P(1), such 

that in the C'^-sense, 

<h£(a;) = $o(a;) + 0{e^), <S>o{x) = ^(Ima;)^ 

and such that the operator 

P, = 0{1) : T-^H^^{\lmx\ < 1/C) T-^H^^{\lmx\ < l/C) (7.4) 

is, microlocally near x {■C = 0})) unitarily equivalent to an operator Pg, 

given in (15.31) . (15.41) . acting on Here 

T : L2(t 2) ^ 77,j,^(CV27rZ2) 

is the standard unitary FBI-Bargmann transform on the 2-torus, associated to the 
quadratic phase function i{x — y)"^/2, as discussed in [101, and we have written 

= }lol{Q)nL\Q,e-^‘^^/’^L{dx)), 

for H C C^/27rZ^ open, including the Floquet periodic versions of the spaces. Let 
us also point out that the unitary equivalence between the operators P^ in fl7.4p and 
Pg is realized by means of a microlocally unitary h-Fourier integral operator in 
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the complex domain, quantizing the canonical transformation in fl7.3p . Similarly to 
fl4.43p . we may illustrate it in a commutative diagram, 


R : T-^Hi 


Ima;| < 1/C) -> T-^H^^{\lmx\ < 1/C) 

Ue Ue 


(7.5) 


Ps ■ Lim 




In particular, according to Proposition 15.11 the leading symbol of is independent 
of X 2 , modulo + e:(R — /(^ 2 ))^)- The subprincipal symbol of is 0{e). 

Remark. From [12], we may recall that writing 

we have = kt o (T*T^), where the canonical transformation kt associated 
to T is given by 

T*T^ 3 {y, v)^{y- iy, y) = (x, 0 e RTl 

Let us also remark that the writing fl7.4p is somewhat informal, and a precise state¬ 
ment is obtained by considering the action of the conjugated operator TP^T~^ on 
the space R<i,^(|Ima;| < 1/C), see also [T2] . 


We obtain a globally dehned IR-manifold A C T*M, which is £-close to T*M ev¬ 
erywhere in the C°°-sense, agrees with that set away from p“^(0), and in a complex 
neighborhood of Aq, it is obtained by replacing 

exp{teHG)oKo\T*T^) 


by 


exp {ieHc) o o Kf\T*T^), 


(7.6) 


which amounts to an C(£^)-deformation A^g iii a neighborhood of Aq. Here we recall 
the holomorphic canonical transformation exp {ieHG), identifying A^g and T*M in a 
neighborhood of Aq, and the real analytic canonical transformation kq in fl4.1l) . given 
by the action-angle coordinates near Aq. The spectral analysis required in order to 
compute the extremal eigenvalues of Pe in Theorem 12.11 will be carried out in the 
globally dehned h-dependent Hilbert space H{A), associated to the IR-manifold A 
by the FBI-Bargmann approach. 
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Recalling Proposition 4.2 and taking into account also Proposition 15.21 eliminating 
the a; 2 -dependence in the normal form by means of a pseudodifferential conjugation, 
we may summarize the discussion so far in the following result. 


Proposition 7.1 There exists a globally defined smooth IR-manifold A C T*M and 
a -Lagrangian torus Aq C A, which is an 0{e)-perturbation of the rational torus 
Aq in the C°°-sense, such that when p E A is away from an -neighborhood of Aq 
in A and 

^ 2 S 

(7.7) 

for C > 0 large enough, then we have 


^ 2 ( 5+1 

lmPfip;h) > emiQ^{Ao) + (7.8) 


Here 0 < S < 1/2 is so small that 3> max(h^/^, The manifold A is 0{e)- 

close to T*M and agrees with it away from a neighborhood of p~^{0) fl T*M. We 
have 


Pe = 0(1) : H{A,m) -E H{A). 


Furthermore, there exists an elliptic h-Fourier integral operator with a complex phase 


U = 0{l) : H{A) ^ 


such that microlocally near Aq, we have 


UPs = {P{xi, hDx, S] h) + R{x, hDx, e; h)) U. 

Here P{xi, hD^, S] h) + R{x, hD,,., e] h) is defined microlocally near ^ = 0 m T*T^, 
the full symbol of P{xi,hDx,e]h) is independent of X 2 , and 

R(x, £i h) = 0(£"+' + (6 - /(&))« + A«+‘), /(O) = 0. (7.9) 


Here N is arbitrarily large but fixed. The leading symbol of P{xi, hD^, e; h) is of the 
form 

piO + O(s^), 

where 


p(0=p(/(6),6) + ^7(0(ei-/(6))', ^?(o)>o, /(o) = o, (7.10) 

and T 9 xi I—)■ q{xi, f{f 2 ),^ 2 ) has a unique minimum, when ^2 £ neigh(0,R), which 
is also non-degenerate. The subprincipal symbol of P{xi, hDx,e] h) is 0{e). 
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Using Proposition 17.11 we shall now discuss a priori estimates for the equation 

{Pe-z)u = v, (7.11) 

when u G H{A,m), v G H{A), and the spectral parameter z G C is conhned to the 
region 

^25 

|Re 2 r| < \mz <em.iQoo{Ao)+ 0{^/eh). (7.12) 

When doing so, following na, m, we shall make use of a suitable partition of 
unity on the manifold A, dehned using Proposition 17.11 and consisting of smooth 
functions satisfying slightly degenerate symbolic estimates. Indeed, the presence of 
such slightly exotic symbols is natural here, as we are dealing with methods based 
on the techniques of normal forms, introducing error terms vanishing to a high order 
along the invariant tori. See also [25]. When quantizing the corresponding symbols 
dehned on A, in the case when M = R^, we follow m and reduce the quantization 
procedure to that of Weyl on the standard phase space T*R^, by means of a 
canonical transformation 

K : neigh(p“^(0),T*R^) —)■ neigh(p“^(0). A), 


such that 

k{X) = X + isHcix) + 

and the corresponding unitary Fourier integral operator with a complex phase map¬ 
ping L^(R^) to H(A). In the case when M is compact, we use the Toeplitz quanti¬ 
zation, following [26] . 

Let us consider a smooth partition of unity on the manifold A, 

l = X + V'i+^2 - (7.13) 

Here x ^ C^{A), = C>(e“^'^™'), m > 0, is a cutoff function supported in an 

neighborhood of Aq intersected with the region where iRePgl < /C. Specihcally, 

we shall obtain y by choosing a suitable function xo ^ C^(T*T^), d'^Xo = 

|a| > 0, depending on ^ only, xo = Ao(0) conjugating the operator xoi^D^) by 
the microlocal inverse of the operator U in Proposition 7.1. In particular, we get, 
using that the subprincipal symbol of Pe=o vanishes, 

= O + O : H{A) -A ff(A), (7,14) 
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The function 0 < e (7°° (A) in fl7.13p satisfies 

m>0, 

and is such that 


-25 


|ReP,(p;/!)| > — 


(7.16) 


near the support of V'l- Finally, 0 < -02 £ 117.13p is such that 

V"^V’2 = m > 0, (7.16) 

and furthermore, 'ip 2 is supported in a region invariant under the iTp-flow, where 


Im Pe (p; h)>e inf Qoo (Aq) + 


-l+2<5 


0 ( 1 )' 


(7.17) 


We also arrange, as we may, so that "02 Poisson commutes with p, the leading symbol 
of Pe=o acting on H{A). 

Let us now return to the equation fl7.1ip . Assume that 6 G (0,1/2) is so small that 

for some hxed p > 0. We can then follow the slightly degenerate parametrix con¬ 
struction for Pg — z, near the support of -i/i, described in detail in Section 4 of m 
and obtain that 

(^( 1 ) 


4’iu \ \ < 


-25 


v\\ + 0{h°°)\\u\\. 


(7.19) 


Here and in what follows the norms are taken in the space H{A). 

When discussing estimates for -^ 2 M, let us notice that ImPe(p; h) = 0(e} on A, and 
near supp'^ 2 , we have, in view of fl7.17p and fl7.12l) . 

^ 1 + 2(5 


Im {Pe{p; h)- z)> 


0{l) 


0{^h). 


Therefore, with a new implicit constant, we get near the support of '02; 

^Im {P,{p;h)-z) > 

provided that the following lower bound on £ holds, 

0il+i5) ^ 


(7.20) 


(7.21) 


(7.22) 
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The lower bound fl7.22p is of the same form as fl7.ip . Using as the natural 

semiclassical parameter and applying the sharp Carding inequality, we get in view 

of (OH), 


-Im ((Pg - z)' 4)2U, i^ 2 u) > 
e 


-2(5 


^-0(l)U|ll(f2«IP-0(A“)ll« 


-25 


> 


0 ( 1 )' 


^lJ 2 u\\^ -Oih°°)\\u\\^, (7.23) 


provided that we strengthen (I7.18p by assuming that 


A<A'. „>0. 


(7.24) 


It follows from (I7.23p that 


. 2 ( S -|-1 


0 ( 1 ) 


'4’2u\\‘^ < 0(1)11 n II \ \'ljj 2 U II + Im {[Pe,^p2]'4’2U,^p2u) + 0(h°°) 11 M 11^. (7.25) 


Here 'ip 2 G C'^(^) has the same properties as -02 and is such that ^’2 = 1 near 
supp (t/’ 2 )- When estimating the commutator [P£,-^ 2 ] in (I7.25p . we get by the Weyl 
calculus, using 07.161) together with the fact that the subprincipal symbol of Pe=o 
vanishes and p and "02 Poisson commute. 


eh 


in, i'2] = [n-o. H + o{^]^o{^] + o{^]^o{^ 




eh 


eh 


(7.26) 


Here we have also used that h^ , in view of (I7.22p . Combining (I7.25p and 

(I7.26p . we get 

V’2n IP < 0(1) 11 n 11 11 ' 1 IJ 2 U 11 + ( ^ ) 11 IP + 0(h°°) 11 u | P, (7.27) 


0 ( 1 ) 
and therefore. 


'02^ IP < 


g(l) | 

riS+2 I 


«’ll" + o(4|llfei<IP + o(ft”)ll«ll", 


(7.28) 


Combining 07.24p . fl7.28p . and a standard iteration argument, we conclude that 

0 ( 1 ) 


i^2U 11 < 


>1+25 


n|| + 0(h°°)||M 


(7.29) 
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Using fl7.13p . fl7.19p . and fl7.29p . we obtain the following a priori estimate for the 
problem fl7.1ip . fl7.12p . 


Il(l-X)«ll<^ll^ll + C?(h~)||n||, (7.30) 

which holds provided that 5 G (0,1/2) and the conditions fl7.22p . fl7.24p are fulhlled. 
In the subsequent analysis, we may therefore concentrate the attention on the region 
supp (x), for the cutoff function y in (I7.13p . 

Let us recall that the function y e (7^(A), V^y = m > 0, in fl7.13p is 

supported in an ^'^-neighborhood of Aq intersected with the region where iReP^I < 
Writing 

{Pe - z)xu = yn + [Pe, y]M, 

and applying the Fourier integral operator P of Proposition 7.1, we get 

{P{xi,hD^,e]h) - z)Uxu = Uxv + U[P^,x]u + Tu. (7.31) 

Here, using fl7.9p . we see that 

T = 0{h^) : H{A) Ll{T^), (7.32) 

where M can be taken as large as we wish, provided that the integer N in Proposition 
7.1 is taken large enough. Furthermore, as discussed above, we may arrange so that 

Py = yoP + 0(h°°):P(A)^P(A), 

where Xo = Xo{hDx,e) is of the form 

Xo{^,e) = Xi (I) Ai 

where Xi ^ C'^(R) is a standard cutoff to a neighborhood of 0. In particular, using 
fl7.10p we see that the support of Xo is contained in the region where 

|e| = (P(£^), \p{0\<O{e^% 

Modifying the operator T in fl7.3ip slightly, we get 
{P{xi,hD^,e-,h) - z) Xi A 2 + U[P^,x]u + Tu (7.33) 
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In the subsequent analysis we shall therefore be working on the cotangent space 
T*T^, in the region where 

6 = (7.34) 

while 

^2 = 0(e^^). (7.35) 

Taking a Fourier series expansion in X 2 , we get a direct sum decomposition 

P{xi,hD^,e;h) = ^2 = h{j - 62 ), (7.36) 

i6Z 

where, according to fl7.35l) . the summation is restricted only to those j & 'Zi for which 
^2 = We shall consider the question of inverting the operator 

F(xi, hD^, e-,h) - z = ^ hD^^,^ 2 , £; h) - z), (7.37) 

j 

where, compared to fl7.12p . the real part of z will be localized further to the region 

|Rez| < (7.38) 

C ye 

where C > 0 is large enough but hxed. Since in Proposition 7.1 we have introduced 
errors that are 0{h^), M 3> 1, see (I7.32p . we would first like to show that the 
one-dimensional non-selfadjoint operator 

P{xi, hDx^^^2i S]h) — z ■. Lg^(T) —)■ Ll^(T) (7.39) 

is invertible, microlocally in the region where .^1 = 0 {e^), with an inverse of tem¬ 
perate growth in 1 /h, when ^2 = is such that |.^ 2 | > h/C\^/e, for a suitable 

hxed Cl, satisfying Q < Ci < C. In doing so, it will be convenient to distinguish 
two cases, depending on the sign of ^ 2 - 

Case 1. Let us assume hrst that ^2 = is such that 

6 > (7.40) 

Then, after a unitary conjugation, we can write, on the level of operators, 
e-*7(«2piFp(a;i, hD^^,^2, S] h)e*7(«2)^iA _ ^ 

= P{f{^2),^2) + g{f{^2) + hD,,,i2){hD,,f + ieq{xi, /( 6 ) + + 0{e^) 

+ hO{e) + 0{h^)-z. (7.41) 
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Here the conjugated operator, acting on the space of Floquet periodic functions 
is still considered microlocally in the region where = 0 {e^), since 
/(C 2 ) = ^>{^ 2 ) = 0{e'^^). Recalling that the derivative of the function ^2 

is strictly positive near ^2 = 0, we conclude, using fl7.38p . fl7.40p . and the 
positivity of g{f{^ 2 ) + ^ 2 )(h-Dxi)^, that the real part of the operator in (I7.4ip . 

which is of the form 

p(/(^ 2 ),^ 2 ) + g{f{^ 2 ) + ^ 2 )(h-Dxi)^ — Rez + 0{e'^ + h^) + hO{e), 

is > h/{C^/F)^ for some C* > 0, and is therefore invertible, microlocally in the region 
^1 = 0{e^), with the norm of the inverse being 0{\/e/h). Here we also use that 
<C h/^/e, in view of fl7.ip . It is therefore clear that the full operator in fl7.4ip 
is invertible, microlocally in the region ,^1 = O(e^), with a microlocal inverse of the 
norm 0 {y/e/h). 

Case 2. We assume now that ^2 = is such that 

Similarly to fl7.4ip . we write 

hD^^,^ 2 , e; - z 

= dif {^ 2 ) + hDxi,^ 2 ){hDxi)^ + ieq{xi, f {^ 2 ) + h-Da:i, >^ 2 ) + C>{e'^) 

+ hO{e) + 0{h^)-w, (7.43) 

where 

^-p(/(6),6) (7.44) 

satishes 

Imtc < £inf Qoo(Ao) + C(\/7h), (7.45) 

(-' 2 V^ 

for a suitable C 2 > 0. In view of (17.351) . we have 

f{^2 ) + 'Cl! ^ 2 ) = 0) + 0{e‘^^), (7.46) 

and therefore, on the operator level we obtain that 

= gif {^ 2 ) + hD^^,^ 2 )ihD^^f + ieq{xi,hD^^,0) + R-w, (7.47) 
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where 


R — + eh + + h^) : (7.48) 

We may also assume that in (I7.47p . the operator q^xi^hD^^^R) is given by the 
classical h-quantization. It follows from fl7.33p that thanks to the presence of the 
cutoff Xi{hDx/s^), to invert the operator in fl7.39p . microlocally in the region where 
^1 = O(e^), we should consider the following equation 

(fi'(/(6) + hD^^,^2){hD^^f + ieq{xi,hD^^) + R-wi) xi u = v, (7.49) 

for u, V E Here Wi = w — ieird Qoo(Ao) and 

6) = 0) + /c(xi, ^i)(p , 

where 

g(xi, 0) = {q) 2 {xi, 0) - inf (5oo(Ao) > 0, 
k{xi,^i)=^i (%^ (xi,f^i) df 

Jo 

and (p G C“(R) is such that ip = 1 near supp(xi)- In particular, the function 
g(xi, ^i) satishes the assumptions for the function V in Proposition A.4 in Appendix. 

Let us set 


H(xi, hD^^) = fi'(/(6) + hD^^,^2){hD^^f + ieq{xi, hD^^), 

We would like to invert the operator A{xi, hD^^) + R — Wi, occurring in the left 
hand side of fl7.49p by an application of Proposition A.4 and to that end, we shall 
assume that 

< h. (7.50) 

Write 

H(xi, hDxi )-wi = e ^ ^ ^ (7.51) 

\ s J e e 

It follows from fl7.45p that the spectral parameter W 2 satishes 

ReW2 > ^ , ImtC2<0(h). (7.52) 
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Here we recall that h = h/^/e. In order to be able to apply Proposition A.4 to fl7.5ip 
we finally have to impose the smallness condition 

(7.53) 

and using (I7.35p . (I7.38p . (I7.44p . and (I7.5ip . we see that (I7.53p holds provided that 

-e^ -C 1. 

£ 

We shall therefore require that the condition 

h < (7.54) 

holds. Once the conditions fl7.5Up and fl7.54p both hold, we are in the position to 
apply Proposition A.4 to fl7.5ip . obtaining that 

{A{xi,hD^^) -wi)~^ 

= s-'0(h-^/^ : i2.+/,6)(T) ^ (7,55) 

We get, using that |ta 2 | > h/{C 2 e^^‘^), 

(A(xi,hD^,) - w,)~‘ 

= e-‘0(h-V'^) = ^ (7,56) 

Returning to the equation (17.491) . we would like to use a standard Neumann series 
argument to invert the operator A{x\, +R — wi in the left hand side of (17.491) . 

and according to (17.561) and (17.481) . we know that this is possible provided that 

(£2-5+1 +£h + e^ + h^) < 1, (7.57) 

which, in view of (17.11) and (I7.22p . is equivalent to the condition 

h-Vi+2'5 < 1. (7.58) 

Comparing the upper bounds (I7.58P and (I7.50p . we see that the latter is implied by 
the former, provided that 

0 < 5 < 1/9. (7.59) 

In what follows, we shall adopt the smallness condition (I7.59p . We arrive therefore 
at the following upper bound on £, 

£«h6/(5+125), (7.60) 
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which is a strengthening of the upper bound in fl7.1l) . 


We shall now also examine the lower bounds on e that we have imposed in the course 
of our argument in this section. To that end, we recall that the lower bounds have 
been introduced in fl7.22p . fl7.24p . and fl7.54p . Comparing hrst the lower bounds 
fl7.22p and fl7.54p . we see that we have 

0(1+4S) ^ j^l/il-5) 

when fl7.59p holds, and our lower bound on e: becomes 




(7.61) 


We should then check the validity of fl7.24p when fl7.6ip holds, and to that end we 
observe that indeed. 



thanks to fl7.59p . 


Combining the bounds fl7.60|) and fl7.6ip . we get the permissible range 

hm-S) 


(7.62) 


where 6 G (0,1/9). The range in fl7.62p is non-empty for 6 G (0,1/9) precisely when 


1 ^ 6 
1-5 ^ 5^ 



1 

18 


Let us summarize the discussion above in the following result. 


Proposition 7.2 Let us consider the operator 


P{xi,hD^,e-,h) = ^P{xi,hD^^,^ 2 ,£-,h), ^2 = h{j - 6 * 2 ) = 


microlocally in the region ^2 = where 1/18 < 6 < 1/9. Assume 

that the spectral parameter z E C is such that 



(7.63) 
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for some constant C > 0. Assume furthermore that 




(7.64) 


and let us assume that the quantum numbers ^2 = satisfy 



Then there exists a family of operators 




such that 


{E{f2,£;h){P{xi,hD^^,f2,£;h) - z) - 1 ) x 


(^) = oun: Li ^ l : 


2 

'^1 


for every x G C'^(R) with support in a sufficiently small but fixed neighborhood of 


0 . 


Remark. Notice that to reach powers of h that are < 1 in fl7.64p . it suffices to take 


5 > 1/12. To obtain the range in fl7.64l) that is as large as possible, we should choose 
5 G (1/18,1/9) to be close to 1/9. 

In what follows, we continue to assume that the spectral parameter G C is conhned 
to the region fl7.63l) . and we shall assume that fl7.64p holds, for some 5 G (1/18,1/9). 
It follows therefore from Proposition 17.21 that in the orthogonal sum decomposition 
fl7.37p . we can restrict the attention to the quantum numbers ^2 = h {.3 — 6 * 2 ) such 
that 



(7.65) 


Using this rehned localization in the parameter f 2 ^ we shall now proceed to show 
that the spectrum of the operator Pg in the region fl7.63|) is contained in the union 
of the pairwise disjoint bands of the form 


b(/( 6 ), 6 ) - Rez| < Co^/eh, ^2 = h{j - 62 ) = 0{h), (7.66) 


where Uq > 1 is large enough but hxed. When doing so, we shall proceed similarly 
to the arguments above, relying upon Proposition A.4 and treating the parameter 
.^2 in a perturbative way. 
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Let us assume that z E C satisfies fl7.63p and is such that for some sufficiently large 
fixed Co > 1, we have 

b(/(6),6) - RezI > Coy/eh, (7.67) 

for all ^2 = h{j — 62 ) = 0{h). Similarly to fl7.43p . we write 

e-R(«2RR^P(a;i, 6, £; - z 

= gU{i2) + hD^,,i2){hD^,f + ieq{x^J\i2) + hP.,,6) + 0{e^) 

+ hO{e) + 0{h:^)-w, (7.68) 

where w = z —p{f{^ 2 ),^ 2 ) satisfies 

|Retc| > Co\/ih, Imw < e: inf (5oo(Ao) + C(\/£h). (7.69) 

Now, in view of fl7.65p . we have 

Q{xi,f{^ 2 ) + hP>^i,6) = q{xi,hD^^, 0 ) + 0 (h), 

and arguing as in the discussion of ’’Case 2” above, we see that we have to invert 
the problem 

(5'(/(6) + hD^^,^ 2 ){hD^J‘^ + ieq{xi,hD^^) + R-wi) Xi {-^P^ u = v, (7.70) 

where q{xi, hD^^) satisfies the assumptions in Proposition A.4 and 

R = 0{eh + + eh + h^) : -t Lg^^y(^ 2 )(T). (7.71) 

The spectral parameter tci in fl7.70p satisfies, in view of fl7.69p . 

- iRewiI > Coh, -Imtci < C>(h). (7.72) 

An application of Proposition A.4 gives, as before, 

(s(/(6) + hDr^,i2)(hD^X + iEq(xi,hDrCl - ""i)”' 

= E-'O k,r‘'’£‘''=') : ii.+/te)(T) ^ (7.73) 

and using fl7.72p . we see that the bound on the operator norm in fl7.73p does not 
exceed 

e-^0{Cp^^h-^). (7.74) 
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To invert the full operator 


5 '(/(6) + + ieq{xi, hD^^) + R-wi, (7.75) 

in the left hand side of fl7.70p . in view of fl7.7ip and fl7.74p . we have to check that 

(^eh + <C 1, (7.76) 


which is satished for Cq > 1 large enough, since clearly, e h, in view of fl7.64p . 
The bound on the norm of the inverse of the operator in fl7.75p is therefore also 
given by fl7.74p . 

Combining Proposition 7.2 with the discussion above, including the estimates fl7.73p . 
fl7.74p . we conclude that if 2 ; G C satishes fl7.63p and is such that fl7.67p holds, 
then the operator P{xi, hD^, e]h) — z is invertible, microlocally in the region where 
^1 = ^2 = C^(^^'^); with a microlocal inverse of the norm 

{P{xi, hD^, £; h) - z)-^ = ^ LI (7.77) 

Coming back to (I7.33p . we therefore obtain for such z’s, 

(7.78) 

where M' 3> 1. Here we have also used fl7.14p . Combining (17.7811 and fl7.3np . we 
obtain the following result. 

Proposition 7.3 Assume that 

< e < /^6/(5+125)^ 


Uxu\ \ < 0{e-^^^h-^)\\v\\ + 


{ sh\ 


for some 6 G (1/18,1/9). Then the spectrum of the operator P^ : H{A,m) —)■ H{A) 
in the region 


Re^;! < 


h 


Imz < einf Qoo(Ao) + 0 {y/eh). 


is contained in the disjoint union of the bands of the form 

b(/(6),6) - Re^l < Coy/eh, ^ = h{j - 62 ) = 0 (h), 


(7.80) 


(7.81) 


where Cq > I is large enough but fixed. 
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We shall finally obtain a precise description of the spectrnm of in the region 
fl7.8ip . for a given value of j G Z, such that ^2 = h{j — 6*2) = 0{h). In doing so, in 
view of the localization for Imz, we may assume that 


l^-p(/( 6 ), 6 )-*^(?) 2 (a^i( 6 ),/( 6 ), 6 )l < Coy/eh, (7.82) 

where ^2 = h{j — 62 ) = 0 {h), and we then know that only the operator 

P (xi, —)■ 

in fl7.36p contributes to the spectrum in this region. Let us introduce the 
elliptic operator 

Q((.A;?2) =9(/(&),6)C? + ^(9^.(«>2(xiK2),/(6),f2))(", (7.84) 

and let ek ,^2 ^ A: G N, be eigenfunctions of Q{t, Dt] ^ 2 ) corresponding to the 

eigenvalues \k{,i 2 ) given in (16.431) . Let also fk ,^2 be eigenfunctions of the adjoint 
Q*(t, Dt, ^ 2 ), corresponding to the eigenvalues Xk{^ 2 )- An application of Proposition 
6.2 allows us to conclude that if fl7.82p holds and the rescaled spectral parameter 

{z - p(/( 6 ), 6 ) - i£{q) 2 M^ 2 ), /( 6 ), 6 )) 

avoids a small but fixed neighborhood of the eigenvalues Aa;(^ 2 ) in the disc \z\ < Cq, 
then z is not in the spectrum of the operator in (17.831) . with 

{P{xi,hD^^,e;h) - z)~^ = O ^ 

In view of the analysis above, we conclude that then ^ Spec{Ps). It remains 
therefore for us to discuss the setup of the global Grushin problem for when the 
spectral parameter z G C is such that 

;^^ 2 ;-p(/(^ 2 ), 6 ) -^£(?) 2 (a^i( 6 ),/( 6 ), 6 )j e neigh(Afc(^ 2 ),C), (7.85) 

for some /c G N with k = 0{1). Using the notation of Proposition 7.1, let us set 

P+: P(A)C, P+M = P+(^ 2 ,fc) (Uyu, 652)^2 , (7.86) 


(7.83) 

quadratic 
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where e^ 2 (^ 2 ) = and R+{^ 2 ,k) has been introdnced in fl6.25l) . nsing the 

eigenfnnctions fk,^ 2 - Dehne also 

R_:C^H{A), R_u-= U-^ {R-{^2,k)u-®e^,). (7.87) 

Here R-{^ 2 ,k) has been introdnced in fl6.25p . and U~^ is a microlocal inverse of 
U. Argning as in Section 6 of [3, we obtain that when fl7.85l) holds, the Grnshin 
operator 

V{z) = : H{A,m) x C ^ H{A) x C 

is invertible, and the corresponding effective Hamiltonian E_^{z) : C ^ C vanishes 
precisely when z is of the form fl6.42l) . fl6.43l) . This completes the proof of Theorem 

2 . 1 . 

8 Numerical illustrations of spectra 

The pnrpose of this section is to present the resnlts of nnmerical compntations of 
the spectra of P^, in the following sitnation, which is easily implemented: let ns 
consider 


Re = + ieq{x, y] hD^, hDy), 

q{x, y; hD^, hDy) = qo{x, y) + qi{x, y)hD^ + q2{x, y)hDy, 

on the torns M = T^ ^ = (R/27rZ)^. Here qo,qi,q 2 are real trigonometric polyno¬ 
mials of degree < F G {1,2,...}. We shall consider the spectrnm of this operator 
near the energy = 1- 

The general assnmptions (12.61) . (12.71) . (12.9p are fnlhlled, the operator P^=q is selfad- 
joint, and the leading semiclassical symbol is of the form (I2.1UI) with 

P = + (l{x,y,^,v) = qo{x,y) + qi{x,y)^ + q2{x,y)r]. (8.2) 

We also have dp ^ 0 along p“^(l) fl T*M. 

The Hamilton flow of p is completely integrable and we have the decomposition 
(12.lip , for p“^(l) rather than p“^(0), where 

'^= U = T^,y X 

(b»?)6T 
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We have 


g!{x,v)= Y, 

\j\,\k\<F 

where the reality of is equivalent to the property, 

Qe{-j,-k) = qi{j,k). 

Rather than taking some particular explicit choice of q, we generate % at random 
by choosing 

k) = ^ {MJ: k) + Ai{-j, -k)) , Ai{j, k) = 

where Oi^j k ^ -^(0)1) ^ire independent Gaussian random variables. The parameter 
K > 0 induces an off-diagonal exponential decay, corresponding to the assumption 
that q{x,y;^,ri) is analytic in {x,y). Then, 

q{x,y,^,v)= 

(i.fc)e[-F,F]2 


where 

qU, k] y) = qo{j, k) + gi(j, k)^ + gaO', k)y. 

Here and below it is understood that [—F, F] is an interval in Z. 

Let G J be a rational torus, so that (^,//) G T and ^/y G Q U {cxo}. The 
Lfp-trajectories in are of the form 

7 ; R 9 s ^ ((xo, yo) + 2s(^, y), (^, y)). 

The restriction of q to such a trajectory is 

qhis);i,y) = Y 1 ( 8 . 3 ) 

(jA)6[-F,F]2 


For the corresponding limit of the trajectory average. 


(q)^ = lim — 
r^ooT 


-r/2 


-T/2 


qh{s);^,y)ds, 


the terms in fl8.3p with y) ■ (j, fc) 7 ^ 0 give a zero contribution, and we get 

(<(>7= E 
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If we write {^, 1 ]) = {—n,m)/\{m,n)\ with {n,m) G Z^, gcd(n,m) = 1, then n 
= 'Zi{m,n) and the intersection of this set with [—F,F]‘^ (viewed as a subset 
of Z^) is equal to 


{/i(m,n); fi E Z, |/i| < F/niax(|m|, |n|)} . 


This gives, 


[F/ max(|m|,|n|)] 


(«)t= F, W™. 


(8.4) 


-[F/max(|m|,|n|)] 


where ^( 7 ) = (a:o, Vo)' {'^1 n) varies in [0, 27r[ and can take any value in that interval 
and [•] denotes the integer part. It follows that 



When max(|m|, |n|) > F this interval reduces to the torus average $'(0,0;.^,p) = 
(g)A^^, so we get non-trivial intervals only for the hnitely many values {m,n) G 
[—F,F]‘^ with gcd(m,? 7 ,) = 1. 

We have written MatLab programs for the production of q and for the calculation 
of (?)a 5 Qoo(A), as well as the supremum and inhmum of q over each torus in J. For 
the graphics, we parametrize J by arg (.^ + ir]) and the hgure below shows: 

• The torus average {q)A, 

• The torus max and min of q, 

• Qoo(A) for each relevant rational torus. 
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q and averages F=2, kappa=2 



arg(xi+i*eta) 


By running the simulation several times we get a series of figures where quite a few 
exhibit the features above. In order to have a numerical illustration of the main 
result of this work it is important that some of the vertical segments (corresponding 
to Qoo(A) for rational tori) reach above the supremum or below the inhmum of the 
curve of torus averages. A larger F will produce a richer picture with more vertical 
segments, but it will also complicate the numerical calculations of the eigenvalues, 
so we settled for F = 2 as a reasonable choice. We also found that k = 2 produces 
some - not too many - visible vertical segments. 

Once an interesting q has been selected, we compute the spectrum numerically 
by working on the level of Fourier coefficients. Thus, if we are interested in the 
eigenvalues with real parts in [Fi, F 2 ], where Fi < F 2 are close to 1, we work with 
Fourier modes for (j, k) in the set 8 of (j, k) G Z^, satisfying + {hkY G 

[Fi,F 2 ], i.e. 

\{j,k)\ G [v^/h, v^/h]. 

The number 4j^8 of such modes is ~ 7r(F2 — Ei)/h? and we ask MatLab to compute 
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the spectrum of the £ x £^-matrix Ae 


given by 


= {ae{j,k;j,k) 


ij,k), {j,k)&£ 


adjAA^k) =h {j + k 


k) 

+ - Ak-k)+ qdi - jA- k)hj + %{j -j,k- k)hk) . 


We cannot let be larger than a few thousand and still we would like h to be 
small and the energy shell thick enough so that the eigenvalues with real 

part inside are not influenced by boundary effects. In the simulations below we 
have chosen the same q as the one in the hgure above and we settled for Ei = 0.85, 
E 2 = I, h = 1/100, leading to ~ 5000. Since the spectra are of width e, we 
rescale the imaginary axis and represent graphically the set of (Re 2 ;, \m.z/e) for ^ in 
the spectrum of P^. We let £ take the values h/2, h, 2h, 4h, 8h, 16h, in agreement 
with Theorem 2.1. 


Spectrum of p + i*epsilon*q,, epsilon=0.005, h=0.01, kappa=2, F=2 



0.5'-^^-' 

0.85 0.9 0.95 1 

Re z 
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Im z/epsilon 


Spectrum of p + i*epsilon*q,, epsilon=0.01, h=0.01, kappa=2, F=2 



0.5'-^^-' 

0.85 0.9 0.95 1 

Re z 
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Im z/epsilon 


Spectrum of p + i*epsilon*q,, epsilon=0.02, h=0.01, kappa=2, F=2 



0.5'-^^-' 

0.85 0.9 0.95 1 

Re z 
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Im z/epsilon 


Spectrum of p + i*epsilon*q,, epsilon=0.04, h=0.01, kappa=2, F=2 



0.5'-^^-' 

0.85 0.9 0.95 1 

Re z 
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Im z/epsilon 


Spectrum of p + i*epsilon*q,, epsilon=0.08, h=0.01, kappa=2, F=2 



0.5'-^^-' 

0.85 0.9 0.95 1 

Re z 
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Spectrum of p + i*epsilon*q,, epsilon=0.16, h=0.01, kappa=2, F=2 



0.5'-^^-' 

0.85 0.9 0.95 1 

Re z 


These eigenvalues form a kind of a centipede with legs sticking out from the main 
body. The majority of the eigenvalues are in the body whose position corresponds 
nicely to the range of the curve of torus averages on the hrst picture. The legs reach 
out to the supremum of the highest and the inhmum of the lowest vertical segments 
corresponding to (5oo(A) for rational tori A. 

Undoing the scaling of Im^, the inclination of the legs should theoretically be close 
to 45 degrees and by measuring this for one of the legs on one of the hgures we 
found some (but not excellent) agreement. 

The main result of this work, Theorem 2.1, describes the individual eigenvalues near 
the extremities of the legs in terms of rational tori. A mathematical treatment of the 
eigenvalues further inside seems more difficult because of the pseudospectral effects 
that are likely to get stronger there. 

By staring at the pictures directly from the pdf hie and creating a movie by switching 
the pages, we see that most of the (rescaled) eigenvalues remain hxed while those 
in the legs and some others move. The hxed ones probably correspond to irrational 
tori and the moving ones to tori that are rational. 
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A Subelliptic estimates for Schrodinger type ope¬ 
rators 

The purpose of this appendix is to establish suitable resolvent estimates for some 
non-selfadjoint operators of Schrodinger type, instrumental in the pseudospectral 
analysis of Section 7. While in the considerations of Section 7, we are concerned 
with operators on the one-dimensional torus, it will be convenient to analyze the 
case of R hrst. See also [T5] . 

Let 

Fo = g{hD^){hD^f + iV{x), V e C'°°(R;R). (A.l) 

Assume that the function g G C'°°(R; R) is such that 


^7-leCo~(R) 

(A.2) 

with 

^>1, ie<7'(oi<i- 

(A.3) 

We mav notice that the conditions llA.3|l are invariant under the sealine offi — >■ 
g{\C), A > 0. We also assume that the potential V is such that 

0<y diVeL°°{R), j>2. 

(A.4) 

and let us make the ellipticity assumption. 


Vix)>^, \x\>C, 

(A.5) 

for some constant C > 0. The semiclassical symbol of Pq, po{x,^) = + iV{x) 

satishes po ^ S{m), where m{x,^) = l+x‘^+^‘^, and when equipped with the domain 
H{m), the natural Sobolev space associated to the order function m, the operator 
Pq becomes closed densely dehned on L^(R). The spectrum of Pq is discrete and we 
notice that 

Spec (Fo) C Po(R-^) = {z e C;arg 2 r G [0,7r/2]}. (A.6) 

Let us make the basic assumption that 


i/-i(0) = {0} C R 

(A.7) 

and that 

l/"(0) > 0. 

(A.8) 
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We are interested in estimates for the resolvent of Po, 

when the spectral parameter z G C is such that 0 < Imz < 0{h) and \z\ 3> h. 
When establishing those, we shall combine some of the results and techniques of |1] 
and [6]. 


In what follows, rather than working with Pq, h will be convenient to consider the 
operator 

P = -^Po = y(a:) - tg{hD,){hD^y (A.9) 

with the symbol p = pi + ip 2 , where 


pi(x,0 = = Vo{x) + Vo(x) = -V"{0)x^ > 0, 


and 


P2{x,^) = -^( 0 ^^- 

It follows from flA.Sp that 

\d^P2\ ~ 1^1 , 

and therefore we obtain the fundamental property 

Vo{x) + HlVo{x,0 ~ I (a:, 01% (^,0 e 

Following [1], let us set, writing X = (a;,0 ^ R^i 


(A.IO) 


(A.11) 


Here M > 1 is a constant to be taken large enough and -0 G ^^(R; [0,1]) is such 
that supp (-0) C (—2,2) and -0 = 1 on [—1,1], It is then straightforward to verify 
that in the region where |W| > we have 


h2/3 

Ga = 0{h), Hc, = 0(l)-^ = 0{h'G), 


and 


d^Go = 0{1) 


1^2/3 f^l/3 


h2/3 


|A|‘''=> |.A| 


H-^ H- -tr 1 ~ ^(1)- 

I W|2/3 I wd/3 ; 


1^1 


(A.12) 


(A.13) 
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Indeed, the validity of flA.12p and flA.lSp follows easily once we observe that in the 
region where 0 < MVq{x) < 2{h and |X| > we have 

H,,Vo = Oi\X\h^/^\X\^^^), XiHp,Vo) = Oi\X\), X\Hp,Vo) = 0(1 + \X\), 


and 


VM ^ M 


Vo^\ 

{h\x\y/^) 


= 0(1) 


hj/3 \xf^ 


, J = l,2. 


Still working in the region |X| > and following [Ij closely, let us obtain a lower 
bound for the function Vq + eoHp^Go, where £o > 0 is a constant to be chosen. We 
have, in view of flA.12p . 

and therefore, in the region where MVq > we get 


Vq + SqHp^Gq > 




M j ^ - 0{l)M 

if we choose £o > 0 small enough. In the region where MVq < we have 

_ h2/3^P2^0 

^0- h ,, 4 / 3 , 


IX 


and therefore 


where 


w + e^Hp,G^ = Vo + Soh^'^ |xr'/' Hl^V^ + R, 


R = I.YI-/’ = O = O (^) 

Using flA.lOp . we see therefore that 


y , .fj r p( ^0 \ ,,2,3 I ,.|2/3 ^ 

Vi + Eoi/^GoSEo—|X| > EO . 

provided that we take M sufficiently large but hxed. It follows that in the entire 
region where |X| > we get 


Vo + EqHp^Gq > 


C^(l) ■ 


(A. 14) 
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We shall now extend the dehnition of Gq to all of R^, and following |1], let us set 

X 


G{X-,h)= il-x 


hl/2 


GoiX-,h). 


(A.15) 


Here x ^ [O 5 1]) is such that x = 1 when |X| < 1. It follows from flA.12p 

and (1A.13P that 

G = 0{h), Hg = d^G = 0(1). 

Furthermore, using flA.ldp we immediately check that on all of R^, we have 

^ 2/3 \Xf^ 


Vo + ^oHp^G > 


0 ( 1 ) 


-0(h). 


We may therefore summarize the discussion above by stating that there exist con¬ 
stants £0 > 0 and Co > 0 , such that we have for all h > 0 sufficiently small, 


c2/3 |jy|2/3 

l/o(:r) + eoHp,G{X) + coh > , X = (x, 0 G R^ 


(A.16) 


Here the real-valued weight function G = G{X,h) G C'°°(R^) has been dehned in 

(ET]i . (IAhsI) . 

It follows from flA.Sp . flA.7p . flA.Sp that pi(X) = V{x) > {1 /G)Vq{x)^ for some 
constant O > 1, and therefore using flA.lbp . we get with ei > 0, ci > 0, 


Rep(X) + e^H,^pG{X) + cih > ^ {h \X\f '^, X e R^. (A.17) 


The estimate flA.17p is analogous to the estimate (4.26) of Section 4 of | 6 ], if we take 
fco = 1 there. Taking flA.17p as the starting point and arguing exactly as in that 
work, we hnd that everything works without any change, provided that the spectral 
parameter 2; G C is such that for some hxed Cq > 1, we have 

Rez < C>(l)h2/3|z|^/^ Gh<\z\<Go. (A.18) 

Here C 3> 1 is a constant large enough and the implicit constant in (IA.18D does not 
depend on G. We therefore obtain the a priori estimate 

0^\z\^/'^\\u\\l2 <0{1)\\{P - z)u\\l^, mg5(R”), (A.19) 
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for z satisfying flA.lSp . 

It therefore remains to discuss the case when z G C is such that 

Rez <0{l)h^/^\z\^^\ k|>Co. 


(A. 20 ) 


Continuing to follow the arguments of Section 4 in [ 6 ] , we obtain from the equation 
(4.34) there that there exist positive constants ci, C 2 such that for z E C satisfying 
(EOl) . we have 


(F — z)u \ \i ,2 11 M I |l 2 + 




ixr 


M, M ) > 


u\\l2. 


(A. 21 ) 


(A. 22 ) 


Here G [0,1]) is a cut-off near 0 such that 

\P{X)-Z\ > y, 

on the support of V 3 (|X|^ / \A)- The spectral parameter 2 ; in flA.20D can be arbitrarily 
large and when estimating the scalar product in the left hand side of flA.2ip . we can 
apply Lemma 8.2 in |1], exactly as it stands, to conclude that 






M, M < 


C^(l) 


(F — z)u 11^2 + C>(l)h|| M 11^2. 


(A.23) 


Indeed, it is easily seen that the proof of Lemma 8.2 of |1] applies in the present 
situation, using the ellipticity property flA.22p and the fact that the symbol p satishes 

\p{X)\ < 0{l) |X|^, X G R^ G F°°(R2), |a| > 2. (A.24) 

Combining flA.2ip and flA.23p . we get 


Z\\u \\\2 < 0{l)—^\ \ (F — z)u 11^2 + 0(1)11 (F — z)u ||i ;^2 II M || 2 , 2 . 
\z\ 

Here we have written Z = brevity. It follows that 

ZWuW^ <0(l){\ + ^\\(P - z)u\\l,, 


|Z, , 

and using also that Z < |^|, we get 

^ 2/3 II M ||l 2 < 0(1)11 (F — z)u \\l 2 . 


(A.25) 


(A.26) 


(A.27) 


We summarize the discussion above in the following result. 
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Theorem A.l Let Pq = g{hDx){hDxY + iV{x) he such that flA.2p . flA.Sp . flA.4p . 
flA.Sp . (IA.7p . flA.Sp hold. There exist constants Cq > 0 and Ci > 0 such that for every 
C > 1 large enough, we have for 

Imz <cih^/^\z\^^\ \z\>Ch, (A.28) 

the following estimate, 

hf^^\z\^^^\\u\\L 2 <ci\\{PQ — z)u\\L 2 , u & H{m). (A.29) 


It follows that when z satishes flA.28p then z is in the resolvent set of Pq and we get 
the resolvent estimate 


(Po - z)-^ = O (h-^/^ ^ L\R). 


(A.30) 


Remark. The discussion above and the result of Theorem A.l extends to the case 
of operators on R"'. 


In what follows, the result of Theorem A.l will only be applied in the case when 
Imz = 0{h). Furthermore, in the considerations in Section 7, we are concerned with 
operators on the one-dimensional torus T = R/27rZ, and our next task is therefore 
to adapt Theorem A.l to this setting. Let us consider therefore 

P = g{hD,){hD,f + iV{x), (A.31) 

where g G C°°(R;R) satishes flA.2p . flA.3p . and let 0 < 17 G C°°(T) be such that 
17“^(0) = {xo}, with V"{xo) > 0. We may write 

P = {hD,,f+ (p{hD^)+iV{x), G C'o°°(R;R). (A.32) 

Let X G [0,1]) be supported in a small neighborhood of Xq, and such that 

X = 1 near xq. On the support of x, the result of Theorem A.l can be applied and 
we conclude that 

h^/^ \z\^^^\\xu\W< 0(1) 11 (P - z)xu I lz .2 + \\u\\l 2 , n G C'°°(T"), (A.33) 


provided that 


z\ > Ch, Imz < 0{h). 


(A.34) 


Using that on the support of 1 — x, fhe potential 17 is bounded from below and 
Imz = 0{h), we see that for all h > 0 small enough, we have 

Im((P - z){l - x)u, (1 - x)m)l 2 > ^1^11 (1 - X)u\\l 2 , (A.35) 
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and therefore, 


II (1 - x)u\\l^ < C|| (P - z){l - x)u\\l^. (A.36) 

Combining the estimates (1A.33P and flA.36p . we see that for all h > 0 small enough, 
Z\\xu\\l2 + II (1 - x)u\\l^ < 0(1)11 (P - z)u\\l2 + 0(1)11 [P,x]^i ||l2- (A.37) 

Here we continue to write Z = and we notice that h Z. We would like 

to estimate the commutator term in the right hand side of flA.37p . and to that end 
we write, using flA.32p . 


0(1) 11 [P, 11 < 0{h) 11 M 11^2 + 0(h) 11 W 11^2 . 


(A.38) 


The first term in the right hand side of flA.38p can be absorbed into the left hand 
side of flA.37p . and we only have to estimate 0(h)|| hu' \ \l 2 . Now 


{^p{hDx)u,u)L 2 + II hu' 11^2 = Re ((P — z)u,u)l 2 + Rez|| m ||^ 2 , 
and therefore, 

II ||l 2 < II (-P - z)u ll^i^ll M ll^i^ + l^^l^'^^ II ||l 2 + 0(1)11 u 11 ^ 2 . (A.39) 

We obtain, combining flA.37j) . flA.38D . and flA.39j) . 


^11 xw IU 2 + II (1 - x)m ||l2 

< 0(1) 11 (P - z)m 11 + 0(h) \z\^/‘^ 11 M I |l 2 + 0{h) \\{P-z)u\ l^i'l I u I I^i^ (A.40) 
so that 

Z\\xu\\l^ + \\{1-x)u\W< 0(1) \\{P-z)u\\l2 + 0{Z^/^) 11 m | ^ 2 . (A.41) 

Assuming that 


we may absorb the second term in the right hand side of flA.4ip into the left hand 
side, obtaining that 


h^/^\z\^''^\\u\\L2 < 0(1)11 {P-z)u\\l2. 

We may summarize the discussion above in the following proposition. 


76 




























Proposition A.2 Let P = g{hDx){hDxY + iV{x) he such that g — 1 G R) 

satisfies g > I, ^ 1- Assume furthermore that 0 < V E C°°(T), with 

P“^(0) = {xq} and V"{xq) > 0. Then for every C > 1 large enough, we have for 

\m.z = 0{h), l^l > Ch, 


satisfying 

h -C 1, 

the following estimate 


/^2/3 II M ||l 2 < C>(1)|| (P — z)u 11^2. 


(A.42) 


(A.43) 


The a priori estimate flA.43p is equivalent to the corresponding estimate for the 
resolvent of P, which provides a resolvent bound in the model case, fundamental 
for the analysis in Section 7. Now the operators that one encounters there are 
somewhat more general than the Schrodinger type operator in flA.3ip . in that the 
potential V (x) should be replaced by a more general h-pseudodifferential operator, 
which furthermore is multiplied by a small coupling constant. We shall now proceed 
to analyze this more general case, essentially by reducing it to the model situation 
treated above. 


Let us hrst consider the following operator on R, 

P,(x, hD,) = g{hD,,){hD,,f + ieV^{x, hD,fi, (A.44) 

where g satishes flA.2p . flA.3p . and following fl7.ip . we assume that 

-C e: < h^^^. (A.45) 

The function V G C°°(R^) in flA.44p is of the form 


V(x,0 = (^) + e (0,1/2), 


(A.46) 


where V is assumed to satisfy flA.4p . flA.Sp . flA.7p . and flA.Sp . and (p G C'^(R"') is a 
standard cutoff function near = 0. The function k is such that 


d^k G L°°(R2), a G N^, k{x, 0) = 0. 


(A.47) 
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We would like to extend Theorem A.l to the operator Pe{x, hD^), and to that end 
we shall simply inspect the arguments above. Writing 


1 

ie 


P^{x, hD^) 


+ V^{x,^rehD^), h 



(A.48) 


we shall view {l/ie)P£ as an h-pseudodifferential operator with the symbol 


=pi+ip2, 


(A.49) 


where 

Pi{x, 0 = yP2{x, 0 = (A.50) 

Writing 

pi(a;,0) = V{x) = Vo{x) + 0{x^), Vq{x) = ^l/"(0)x^ 

we see that uniformly in e, we have 

Vo + Hl^Vo^\X\\ XeR2. (A.51) 

Here we also notice that d°‘p E L°°(R^) for all a E with |a| > 2, uniformly in 
£. Arguing as in flA.lip . (1A.16D . (lA.17p . we conclude that there exists a real-valued 
weight function G E C'°°(R^) with 

G{X) = 0(h), 


such that for some constants 5i > 0 and Ci > 0, we have for 0 < h small enough, 

Rep(x,0) + 5ihri^pG(X) + cih> ^h2/3|^|2/3^ X = (a:,0eR^ (A.52) 

(_y 

Using flA.46p and flA.52p . we obtain that 

Rep(A) + SiHi^pGi^) + cji > - ly^|g|<o(,.)0(v^ 1^1)- (A.53) 

Here we have used that 


k{x,i) = 0{\i\). (A.54) 

When estimating the right hand side in (lA.53p . we notice that 

^^ 2/3 |^| 2/3 _ ^ - ^V-emOie^)0{y/i\^\^^^)) 

> ^ > 0, (A.55) 
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provided that 


^5/3+l/2-l/6 ^ ^2/3 

The latter condition is equivalent to 


/,2/3 

iw- 


£ 


1+1 


■C h. 


(A.56) 


(A.57) 


Assuming that flA.57p holds, we conclude that 


Rep{X) + 6^Hj^pG{X) + c^h > , X 


(A.58) 


In order to be able to apply the general arguments of [6] and [1] to the operator 
Pe/ ie, similarly to the discussion leading to Theorem A.l above, we should also 
observe that for |^| 3> h, we have in the region where 


|Af < a 

that the ellipticity condition 

\p{X)-z\ > 

holds, uniformly in £. Indeed, this follows from (IA.46p . flA.49p . flA.SOp . flA.54p . 
as well as the fact that £ -C h <C |^|. It is then straightforward to check that 
the arguments in the beginning of the appendix apply and we obtain the following 
result. 


Theorem A.3 Let = g{hDx){hDxY + ieV'^{x,hDx), where -C £ -C 
Assume that g G C'°°(R; R) satisfies flA.2p . flA.3p . V is of the form flA.46p . flA.47p . 
and that flA.4p . flA.SD . flA.7D . flA.SD hold. Assume that 


5 

''2 


< h. 


Then we have, writing h = h/y/e, 


P. 


-1 


= O (h-^P : h2(R) ^ L\R), 


(A.59) 


provided that \z\ > Ch for C > 1 sufficiently large, and 

Im 2 ; < 0{h). 
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Repeating the arguments leading to Proposition A.2, with Theorem A.l replaced by 
Theorem A.3 and with an estimate of the form flA.3511 obtained by an application of 
Garding’s inequality, we next obtain an adaptation of Theorem A.3 to the setting 
of the torus. 

Proposition A. 4 Let Pg = g{hDx){hDxY + ieV^{x, hD^), where <C £ -C , 
and g G G°°(R; R) is such that g — 1 ^ C'^(R) with 


9>h ie<7'(oi<i- 


Assume that V is of the form 




V(x, 0 = V(x) + k(x, f)(f J , 5 e (0,1/2). 

Here 0<Ve C°°(T), V-\0) = {xq}, V"(xo) > 0, k e S(T*T, 1), k(x, 0) = 0, and 
ip e G^(R). Assume that -c h and let us set 

r h 

h = —= -C 1. 

Let z ^ C be such that Imz = 0{h), \z\ > Ch, for C > 1 sufficiently large, and 

-C 1. 


Then we have 


P. 


— z 


-1 


= o( h-2/3 : L 2 ^T) ^ L'^iT 


(A.60) 
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